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Abstract 

Line integration of generalized functions is studied. Second or- 

■ der partial differential equations with piecewise continuous and gen- 

^3 eralized variable coefficients over Cayley-Dickson algebras are investi- 

ctf . 

gated. Formulas for integrations of such equations are deduced. For 
this purpose a non-commutative line integration is used. Examples of 



J> ■ solutions of partial differential equations are given. 

in 



1 Introduction. 

It is well-known, that differential equations have many-sided applications in 
different sciences including physics, mechanics, other natural sciences, tech- 
niques, economics, etc. The differential equations also are very important for 



mathematics [3 [T3J [MJ [T51 123 EZ]- Predominantly differential equations are 
considered over fields such as real, complex, or with non-archimedean norms. 
Recently they are also begun to be studied over Clifford algebras [9| fTU| fTTj. 

Such algebras have a long history, because quaternions were first intro- 
duced by W.R. Hamilton in 1843. He had planned to use them for problems 
of mechanics and mathematics [121 E3] • Their generalization known as the 
octonion algebra was introduced by J.T. Graves and A. Cayley in 1843-45. 
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Then Dickson had investigated more general algebras known now as the 
Cayley-Dickson algebras [U [21 [TJ] . 

The Cayley-Dickson algebras, particularly, octonions and quaternions are 
widely used in physics, but mainly algebraically. Already Maxwell had uti- 
lized quaternions to derive his equations of electrodynamics, but then he had 
rewritten them in real coordinates. 

In the 50-th of the 20-th century Yang and Mills had used them in quan- 
tum field theory, but theory of functions over octonions and quaternions in 
their times was not sufficiently developed to satisfy their needs. Discussing 
that situation they have formulated the problem of developing analysis over 
octonions and quaternions [8]. This is natural, because quantum fields are 
frequently non-abelian |35j. Dirac had used complexified quaternions to solve 
the Klein-Gordon hyperbolic differential equation with constant coefficients. 

This work continues previous articles of the author. In those articles 
(super)-differentiable functions of Cayley-Dickson variables and their non- 
commutative line integrals were investigated [221 E2J Ell ESI 125] • In the 
papers [2H EU] differential equations and their systems over octonions and 
quaternions were studied. 

The Cayley-Dickson algebras A r have the even generator i Q = 1 and the 
purely imaginary odd generators i\, ^-i, 2 < r, i\ = — 1 and ioik = h 
and ikii = —irik for each 1 < k ^ I. For 3 < r the multiplication of these 
generators is generally non- associative, so they form not a group, but a non- 
commutative quasi-group with the property of alternativity ik(ikk) = {"vfyh 
and {iiik)ik = k(i>k) instead of associativity. Ordinary super- analysis operates 
with graded algebras over Abelian groups. Therefore, super-analysis over 
the Cayley-Dickson algebras is in some respect more complicated than usual 
super- analysis, for example, over the Grassman algebras. 

The aim of this paper is in developing of Dirack's approach on partial 
differential equations with variable piecewise continuous or generalized coef- 
ficients. 

The technique presented there is developed here below for solutions of 
partial differential equations of the second order of arbitrary signatures and 
with variable coefficients which may also be piecewise continuous or gener- 
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alized functions. Moreover, signatures may change piecewise in a domain. 
Formulas for integrations of such equations are deduced. For this purpose a 
non-commutative line integration of generalized functions is developed. Ex- 
amples of partial differential equations are given. Moreover the approach of 
§§2-25 over the Cayley-Dickson algebras A v gives the fundamental solution 
of any first and second order linear partial differential equation with variable 
z-differentiable Av- valued coefficients, z G U C A v , where U is a domain in 
A v satisfying some mild convexity conditions described below. These results 
can be used for solutions of concrete partial differential equations or their sys- 
tems of different orders with piecewise continuous or generalized coefficients, 
for example, of Helmholtz' or Klein-Gordon's types, which are important in 
optics of composite materials or quantum field theory. Finally solutions of 
some types of non-linear partial differential equations over Cayley-Dickson 
algebras are studied. 

Main results of this paper are obtained for the first time. 

2 Partial differential equations of the second 
order. 

1. Remarks and notations. For a subset U in either the quaternion skew 
field H = A% or in the octonion algebra O = .A3 or the Cayley-Dickson 
algebra A r , r > 4, we put ir S)P)t (U) := {u : z G U, z = J2 v eh w v y i u = 
w s s + w p p} for each s^peb, where t := Sveb\{s, P } w vV G Ar )S ,p '■— {z G A r : 
z = Eveb^vV, u> s = Wp = 0, Wy G R Vv G b}, where b := {z , n, ^-1} 
is the family of standard generators of the algebra A r so that i| = — 1, 
for each j > 1, ijik = —ikij for each j ^ k > 1, zq = 1. Geometrically 
the domain 7r SjP)t (C/) means the projection on the complex plane C SiP of the 
intersection U with the plane 7f S)P)t 3 t, C SjP := {as + bp : a, b G R}, since 
sp* G b := b \ {1}. Recall that in §§2.5-7 [22] for each continuous function 
/ : U — > Ar it was defined the operator / by each variable z G A T . If a 
function / is z-differentiable by the Cayley-Dickson variable z G U C A r , 
2 < r, then f(z) = dg(z)/dz, where (dg(z)/dz).l = f(z). 
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A Hausdorff topological space X is said to be n-connected for n > if 
each continuous map / : S — >■ X from the k- dimensional real unit sphere 
into X has a continuous extension over R k+1 for each k < n (see also |36j). 
A 1-connected space is also said to be simply connected. 

It is supposed further, that a domain U in A™ has the property that 

(Dl) each projection Pj(U) =: Uj is (2 r — l)-connected; 

(D2) ir s ,p,t(Uj) is simply connected in C for each k = 0, 1, 2 r_1 , s = 12k, 
p = Z2fc+i, t G *4. r ,s, P an d w G C SjP , for which there exists z = u + 1 G Uj, 
where = (0, 0, 1, 0, 0) G .A™ is the vector with 1 on the j-th place, 

Pj(z) = iz for each z G A™, z = YJ£=\ ^ ze ji ~* z e f° r eacri J = 1) •••> m i 
m G N := {1, 2, 3, ...}. Frequently we take m = 1. Henceforward, we consider 
a domain [/ satisfying Conditions (Dl, D2) if any other is not outlined. 

The family of all A r locally analytic functions f(x) on U with values in 
A r is denoted by T-L(U,A r )- It is supposed that a locally analytic function 
f(x) is written in the ^-representation v(x), also denoted by v = . The 
latter is equivalent to the super-differentiability of / (see [231 121] )• Each 
such / is supposed to be specified by its phrase v. 

For each super-differentiable function f(x) its non- commutative line in- 
tegral J 7 f(x)dx in U is defined along a rectifiable path 7 in U. It is the 
integral of a differential form f(x).dx, where 

(11) f(x) = dg(x)/dx, 

(12) [dg(x)/dx)].l = f(x) for each x G U. 

A branch of the non- commutative line integral can be specified with the 
help of either the left or right algorithm (see [231 E21 El]). We take further 
for definiteness the left algorithm if something another will not described. 
For / G T-i(U, Ar) and a rectifiable path 7 : [a, b] — > A r the integral / f(x)dx 
depends only on an initial a = 7(a) and final (3 = 7(6) points due to the non- 
commutative analog of the homotopy theorem in U, where a < b G R. When 
initial and final points or a path are not marked we denote the operation of 
the non-commutative line integration in the domain U simply by / f(x)dx 
analogously to the indefinite integral. 

To rewrite a function from real variables Zj in the ^-representation the 
following identities are used: 
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(1) z 3 = (-zij + ij(2 r - 2)-\-z + ELl 1 ik(zi* k )})/2 
for each j = 1,2,... , 2 r - 1, 

(2) z = (z + (2 r - 2)- 1 {-z + i k (zi* k )})/2, 

k=i 

where 2 < r G N, z is a Cayley-Dickson number decomposed as 

(3) z = zoio+---+Z2 r -ii2 r -i £ »4r, G Rfor each j, i* k = i k = —ik for each 
fc > 0, i = 1, since i fc (i ifc) = i = 1, i*(ijifc) = -h(i* k ij) = -(ki* k )ij = -ij 
for each k > 1 and j > 1 with k ^ j (shortly k ^ j > 1), ik(iki*k) = f° r 
each > 0. 

As usually C°(C/, A,) denotes the R-linear space of all continuous A v - 
valued functions / : U — >■ A v . More generally C n (U, A v ) denotes the R-linear 
space of all n times continuously differentiable by real variables zo, ...,z 2 v-i 
functions /:[/—>■ A v , where n G N. Certainly, C n (U, A v ) can be supplied 
with the structure of left- and right-module over the Cayley-Dickson alge- 
bra A v using point-wise multiplication of functions f(z) on Cayley-Dickson 
numbers from the left and the right. 

2. Factorization and integration of equations. 

We consider the second order partial differential equation: 

(1) Af = g, where 

k k 

A = ai^/drtdTrn + ^aid/dn 

l,m=l 1=1 

is a partial differential operator of the second order. Let us suppose that the 
quadratic form 

l,m 

is non-degenerate and is not always negative, because otherwise we can con- 
sider —A. Moreover, let a matrix of coefficients be real and symmetric 

a 2,m( r ) = a m,«( r ) e R> a h r i £ R- f° r eacn h m — 1, Then we re- 

duce this form a(r) by an invertible R linear operator C = C(r) to the sum 
of squares. This means, that 

fc A; 

(2) A = 5>dW + £/W0*i, 
1=1 i=i 
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where dsj/dri = Cij(t), C = (Cij), with real-valued functions b; and (3[ for 
each /. Here 

bz^M = a p,m,C Pt jC mi i and 

k 

flj = ^ ap >m (dC P j / dT m ) + 2^ a vC v j 

v=l 

for all j, I = l,...,k. In the case when coefficients of A are constant, using 
a multiplier of the type exp(X)z ejsj) it is possible to reduce this equation to 
the case so that if b; ^ 0, then (3i = (see §3, Chapter 4 in [S]). Therefore, 
one can as usually simplify the operator with the help of such change of 
coordinates and consider that only fli may be non-zero if bi = 0. 

Thus one can choose an invertible real matrix (ch, m )h,m=i,...,k correspond- 
ing to C = C(t) so that b ; < for p + 1 < I < k and b ; > for < I < p, 
where < p < k, q := k — p. When q = and 0i = for each / the 
operator is elliptic, for q = and /3i ^ the operator is parabolic, for 

< p < k and /3i = for each I the operator is hyperbolic. Sometimes 
the matrix C can be chosen constant on a domain, where the signature 
(p, q) of the quadratic form o(r) is constant. We suppose that the sums 
Z)f=ibf(x) > and Sf= p +ibf(x) > are positive A-almost everywhere on 
a domain U, where A is the measure induced by the Lebesgue measure on 
the real shadow of the Cayley-Dickson algebra. Generally the natural num- 
ber k — p = q(x) may either be constant or change while crossing the sur- 
face {x E U : Sz=ib; 2 (x) = 0}, when the domain U satisfies Conditions 
1(D1,D2). 

We consider elliptic and hyperbolic partial differential operators reduced 
to the sum of squares 

(3) A=\j2h l (x)d 2 /dx% 

1=0 

where bj(x) G R for all x = xqIq + ... + a^-i^V-i i n the open domain 
U C A r satisfying Conditions 1(D1,D2) in the Cayley-Dickson algebra A r , 

1 < k < 2 r — 1, 2 < r < 3. Practically the coefficient b; can depend only 
on xo,...,Xk remaining z-differentiable in definite ^-representations due to 
Formulas 1(1 — 3) for each 

More generally we can consider partial differential operators of the form 
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(4) A = C\B X + ... + c m B m , where CjBjf = Cj(Bjf), while each 



mi+...+mj 



(4') B 3 = £ Ms^/fti* 



fc=mi+...+mj_i+l 



is an elliptic partial differential operator of the second order by variables 
V+...+m 3 -i+i.->V+...+m J i Cj E A r with Re(cj) > for each 1 < j < I, 
Re(cj) < for every j > I, with \cj\ = 1 for each j = 1, ...m, where 1 < r, 
1 < I < m, m = 0. 

We remind, that Dirac had used complexified bi-quaternions to solve 
Klein-Gordon's hyperbolic partial differential equation with constant coeffi- 
cients appearing in spin problems. That is, he had decomposed d'Alembert's 
operator <9 2 /<9t 2 — V 2 as the product \*a\a over the complexified bi-quaternion 
algebra H c with the first order differential operator a. 

If follow this approach one takes the complexified Cayley-Dickson algebra 

(5) (A r )c =A r ® Arh 

where i is taken to be commuting with ij for each j = 0, T — 1. Now the 
algebra (A r )c is already not the division algebra even for 2 < r < 3, that 
is two non zero elements with zero product occur in it. Then each element 



Let each coefficient Cj be written in the polar form 

(6) Cj = exp(i K(j - ) 7 i ) 
with < |tj| < 7r, j = 1, ...,m, 1 < r, 1 < < n(j + 1) for each j. Put 
p = pi + .--+p m , where pj = for either jj = or = k{j — 1), while pj = 1 
for 7j 7^ and ^ n(j — 1). Up to an isomorphism we take the Cayley- 
Dickson algebra A v with v > r satisfying inequalities T~ l < 2 p (m + 1) < 2 V . 
Further we make the complexification (Ay)c of the Cayley-Dickson algebra 



Take two non-negative integer numbers < r and v with r < v G Z. We 
consider the quotient algebra over the real field A v /A r =: A r>v . For r = v 
this algebra is isomorphic with the real field R. For r < v the algebra A r , v is 
isomorphic with ©^Lo" 1 R^fc- The latter algebra is produced by generators 
{i-2 r i3 '■ P = 2 7 — 1; 7 = 0, 1, ...,v — r} and their finite ordered products, 




-V ■ 
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that gives the generators set {i 2 rk : k = 0,...,2 v ~ r — 1}, where generators 
satisfying the numbering rule iji 2 s = ij +2 s for each 1 < s, j = 0, 1, 2 s — 1 
can be taken up to an isomorphism of the Cayley-Dickson algebra A s +i- 
Therefore, the algebra A r , v is isomorphic with the Cayley-Dickson algebra 
A v - r , since the doubling procedure can be started from another suitable 
purely imaginary Cayley-Dickson numbers such as generators [H [TJ]. But 
we consider in A r , v its specific generators basis {i 2 r k : k = 0, 2 v ~ r — 1}. 

For each Cayley-Dickson numbers x,y G A r we define the real- valued 
scalar product 

(RS) (x,y) = (x,y) r := Re(xy), 
where z = z* denotes the conjugated number, while Re(y) := (y + y*)/2 
denotes the real part of y. 

The real scalar product .) r in A r we extend on the algebra A r , v as 

(SP) <x,y> r:V =xy= ^jyvkivji\r k 

j,k=0 

for each x, y G A TyV , x = YsjLo' 1 ETrflvj-, X2 r j G R for each j = 0, 2 v ~ r — 1. 
Particularly, one gets < x, y >o, v —< x , V >v I n the case of the complexified 
algebra (Ar,v)c the scalar product is: 

(SPC) < (a, b), (c, d) > r ,v=< (a, b), (c, d) >= (< a, c> - < b, d >, < a, d > + < b, c >), 

for all (a,b) and (c,d) G (A riV )c- 

We recall the doubling procedure for the Cayley-Dickson algebra A r+ \ 
from A r - Each Cayley-Dickson number z G A r +i is written in the form 
z = £ + T]l, where l 2 = —1, 1 ^ A r , £, rj G A r . The addition of such numbers 
is componentwise. The conjugate of any Cayley-Dickson number z is given 
by the formula: 

(Ml) z* := C - r]l 
The multiplication in A r+ \ is defined by the following equation: 

(M2) (C + 77l)( 7 + 5\) = (£7 - fy) + (5£ + 777)1 
for each £, 77, 7, 5 G -2 := £ + ??1 G A+i, C := 7 + <51 G A+i- 

Using Formula (M2) we get: (bi 2 r k ) (i 2 r k b)* = {bi 2 r k )(b*i* 2rk ) = b 2 = 
(b 2 i 2 rk)i 2 r k = i2 r k(i2 r kb 2 )* for each k > 1 and 6 G *4. r , since = — ij for each 
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j > 1. Another useful identity is the following: (i s i2 r j)i2 r k = ^i^s^k)^] f° r 
each < s < T — 1 and k ^ j with k > 1 and j > 1, since {i s i2 r j)i2 r k = 
(i s ilr k )i2rj- Certainly also the equality (i s io)i* + (i s ij)i*Q = holds for each 
j > 1 and 1 < s < 2 r — 1, since i$ = 1. Therefore, Formulas (ST 3 ) and (4') 
together with the latter identities imply: 

2 v-r_± 

(6) <cBy,y> v = J2 c < B V^j, Vvj >r,v 

3=0 

for each c <E A r and a twice differentiable function y with values in „4 ri „. 

Relative to the complex scalar product given by Equality (SPC) we de- 
compose the operator A (see (4, 4') above) in the form 

(7) A = (i<r)(i(7i) + Q = -<7<7i + Q, 
where a, o\ and Q are partial differential operators of the first order, each 
complex number a G C is presented as a real 2x2 matrix. Particularly, 
i = g), i* = (° Each subalgebra g Jife constructed from two generators 
ij 7^ i k is associative, consequently, (wi k )(w*il) = w 2 and w{{wi k )i k ) = —w 2 
for each u> = wq + Wjij with w 0? Wj G R» Therefore, we can take 

m m,i+...+m,j 

(8) = ^ 51 a*(^)(9//9^fc)[«'j*i2rfc] and 

j=l k=mi+...+rrij-i+l 
m m 1 +...+rn :j 

(9) <ri/(2;) = X; J2 a k (z)(df/dz 2 r k )[w j i*r k \ 

j=l k=mi+...+mj—i+l 

on the space of ^4 r -valued (super-) differentiable functions / for r < 2 or 
real-valued functions / for 3 < r of the A r , v variable, since (iik) 2 = i 2 il = 
(— l) 2 = 1 for each k > 1, where w 2 = Cj for all j and a|(x) = b k (x) for 
each /c and x, Wj G *4. r , a. k (x) G R for all k and x, i2 r fc G A r>v , z 2 r k = x k , 
z = Efe z 2 r k i 2 r k G df(z)/dz 2 r k = (df ' (z) / ' dz) .i 2 r k ■ For b = d 2 f/dz 2 r k and 
1 = i 2rk and w G A r one has the identities: (b(w\))(w*Y) = ((wb)\)(w*\) = 
—w(wb) = —w 2 b and (((b\)w*)\)w = (((bw)\)\)w = —{bw)w = —bw 2 in the 
considered here cases. The operator Q is given by the equality: 

m mi+...+mjj mi+...+mj 2 

(io) Qf(z)= EE E 

ji J2=l fci=mi+...+m ;)1 _i+l fe2=mi + ...+m 3 - 2 _i + l 

a fcl (*) { (<9a fc2 (*) / dz 2 r kl ) (9// dz 2 . fe2 ) [w h i* 2rk2 \} [w* h i 2 r kl ] , 
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since ii^ = if.i in the complexified Cayley-Dickson algebra (A v )c for each k. 
The latter equality (10) shows, that the differential operator Q is non-zero, 
when a.k(z) are non-constant coefficients. 

If use io = 1 and d/dzo also one can write out d'Alembert's operator in 
our notation 

(11) d 2 /dzl-Y? j=l d 2 /dz 2 j = (i*d/dz +i 1 d/dz 1 +i 2 d/dz 2 +i 3 d/dz 3 )(id/dza+ 

iid/dzi + i 2 d/dz 2 + i 3 d/dz 3 ). 

We recall, that the Cayley-Dickson algebra A r is power associative, that 
is z k z l = z k+l for all natural numbers k and I. But the complexified Cayley- 
Dickson algebra (A r )c is not power associative for r > 3, since the Cayley- 
Dickson algebra A r is not associative for r > 3. Therefore, we do not widely 
use the complexified Cayley-Dickson algebras, but we utilize the Cayley- 
Dickson algebras A v over the real field R, when something other will not be 
specified. 

With these decomposition of operators given by Equations (7 — 9, 11) 
the differential equation (1) can be integrated with the help of the non- 
commutative line integration. We consider at first the partial differential 
equation 

(12) Tf = g 

on an open domain U in A v , where 

(13) Tf^^idf/dz^iz)}, 

f and g and ipj{z) are ^4„- valued functions on the domain U satisfying Con- 
ditions 1(-D1, D2), where g, ipj G C°(U, A v ) for each j, particularly they may 
be A v (super-)differentiable functions. 

3. Line integration over Cayley-Dickson algebras. Take any phrase 

(1) K Z ) = Em{c m ,Z m } q(m) 

corresponding to the function /, where 

{Cm; Z }<j(m) {Cl,m\Z 1 ...Ck^m^Z k "\q(m)i 

q(m) is a vector indicating on an order of multiplications in the curled brack- 
ets, Cj jmj G A v for each j, m = (m 1; ...,m k ), k G N, < rrij G Z for each j, 
z k = (...((zz)z)...)z. We put for convenience z° = 1 in the considered phrases. 
Though the symbol z° can be retained when necessary to specify a branch of 
the line integral over the Cayley-Dickson algebra A r (see [221 ESI [21] ) • Using 
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the shift z i — y i^z — z) we can consider such series with the center at a point 
z instead of zero. Then the derivative of the phrase is: 
(2) dfi{z)/dz = 

^2m,l,j{Cl,miZ 1 ■■■Cj — l,rrij-iZ 1 Cj,rrij ( )Cjj^\^ m . +1 Z i +1 . . .Cfc i?Ttfc 2 fe }<?(m,)5 

where J denotes the unit operator, so that d^/dz is the operator valued 
derivative function, < I < m 3 ■, — 1, j = 1, k. From Equality (2) it follows 
that 

(3) (dfi^px) I dx) .1 = (dfi(z) / dz) .i p = dfi(z)/dz p for z = i p a;. 
If 7 : [a, 6] — > A r is a function, then 

K 6 7 := sup P | 7 (t,+i) - 7(^)1 
is called the variation of 7 on the segment [a, b] C R, where the supremum 
is taken by all finite partitions P of the segment [a, b], P — {t — a < ti < 
... < t n = b}, n G N. A continuous function 7 : [a, b] — > A r with the finite 
variation < 00 is called a rectifiable path. It is convenient to take the 
unit segment [a, b] = [0, 1] using a suitable reparametrization. 

We say that a function v on U is absolutely continuous on U if for each 
rectifiable path 7 : [0, 1] — > A v for each e > and each r e [0, 1] a positive 
number 5 > exists so that V r min ( 1 ' r+<5 )z/( 7 ) < e and V^ ax(0 r _ (5) z/( 7 ) < e. 

We call a function z/ of bounded variation on U if for each rectifiable path 
7 : [a, b] — > U the variation V^u(j) < 00 is finite. The family of all functions 
v : U — > Av of bounded variation will be denoted by V(U, Av). 

The non-commutative line integral J f(z)dv(z) along a rectifiable path 
7 : [0, 1] — > U C A v for a phrase ji and a given function v of bounded 
variation is the limit by partitions P = {0 = tq < T\ < ... < r„ = 1} with 
their diameter S(P) = sup^- |tj +1 — Tj\ tending to zero of integral sums 

J 7 f{z)dis(z) := lim 5{P) ^ j: j (dK(z)/dz)\ z=j{Tj) .[is{~f(T j+1 )) - 1/(7(7}))], 
where (dn(z)/dz).l = fx(z) for all z £ U. The notation 

/(z) = dg{z)/dz and /t(z) = dn{z)/dz 
is also used, where (?(z) is a super-differentiable function to which the phrase 
k corresponds. 

If / is a continuous function we fix for it a sequence f n (z) of super- 
differentiable functions and their phrases fi n (z) such that f n (z) converges 
to f(z) on each compact subset of the domain U, where n G N. The non- 
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commutative line integral has a continuous extension on the R-linear space, 
left and right A v module, of continuous functions C°(U, A v ) for a marked 
function u(z) of bounded variation and a given rectifiable path 7: 
/ 7 f{z)dv{z) := lim^oo / 7 f n {z)du(z). 

This means that the R homogeneous A v additive operator f(z) is defined 
for the continuous function / in the sense of distributions: 

(/; iy,i):=I J f(z)du(z) 

for each rectifiable path 7 in U and every function u(z) of bounded variation. 
Particularly, u(z) = id(z) = z on U can also be taken. 

If v and / are super-differentiable functions such that the derivative 
dv{z)jdz of v is the invertible R homogeneous A v additive operator for each 
z G U, then a super-differentiable solution of the differential equation 

(dg(z) I dz) . (dv(z) / dz) = dq((u(z)) / dz 
on U exists, since dz/dv = (du/dz)^ 1 . That is, (dg(z)/dz).di > (z) = (dq(v(z))/dz).dz. 
Therefore, 

f 1 f(z)du(z) = \\ui s{P) ^Y, j {dq(p)/dp)\ v=v{l{T])) \u( 1 (T j+1 )) - 1/(7(7}))] = 

L(~,)P(y)dy, 

where p(u) = {dq{v)/dv).l (see also Theorems 2.11 and 2.13 in [27]). 

A function v : U — > A v is called piecewise continuous or differentiable or 
super-differentiable on a domain U in the Cayley-Dickson algebra U, if a fam- 
ily of open or canonical closed subsets Uj of U exists so that each restriction 
v\u is continuous or differentiable or super-differentiable respectively, where 
U = \Jj Uj and Uj n U k = dUj n dU k for each j ^ k, dUj := cl(Uj) \ Int(Uj), 
cl(Uj) denotes the closure of Uj in A v and Int(Uj) denotes the interior of Uj 
in A v . 

If / is a continuous function and v is a function of bounded variation for 
which the limits lim n f n = f and v = lim n v n uniformly converge on each 
compact subset of U and phrases fi n of f n and e n of v n are specified, where 
f n are super-differentiable functions and v n are piecewise super-differentiable 
functions on U so that lim„ ^(^"(7) — 2/(7)) = for each rectifiable path 7 
in U, then 

! 1 f{z)dv{z) = lim n ^j n {z)dv n {z) = \im n J vn(ry) p n (y)dy = f v(l) p{y)dy, 
where p(y) = lim n p n (y). This means that under rather general conditions 
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the line integral of the type f f(z)dv(z) relative to the function v of bounded 
variation reduces to the usual non-commutative line integral f v p(y)dy, where 
rj = 1/(7). 

Take the branch of the non-commutative line integral prescribed by the 
left algorithm (see §2 in [2"3"j I2"2"]). The real algebra gk,i, s formed from the 
generators ij, ik and i s is alternative. Each rectifiable path can be presented 
as the limit of rectifiable paths consisting of joined segments parallel to the 
straight lines ijH with respective j. We certainly have {i q ip)i p = —i q for each 
p > 1 and (i q io)io = i q for each q > 0. 

For each j = 0, 2 r — 1 the R- linear projection operator ttj : A r — > Hij 
exists due to Formulas 1(1 — 3) so that TTj(z) = ijZj = Zjif 

(PI) n^z) = (-i^) - (T - 2)- l {-z + J2L- 1 ^k(ztl)})/2 
for each j = 1, 2, T — 1, 

(P2) Tro(^) = (z + (T - 2)- 1 {-z + ^\ fc (^)})/2, 

fc=i 

where 2 < r G N. 

4. Line anti-derivatives over Cayley-Dickson algebras. 

Theorem. Let a first order partial differential operator T be given by 
Equation 2(13) with real-valued continuous functions ipj(z) G C°(U,A V ) for 
each j such that ipj(z) 7^ for each z G U and each j = 0,...,n, where 
a domain U satisfies Conditions 1(D1, D2) , qz is a marked point in U, 
1 < n < 2 V , 2 < v. Then a line integrally : C°(U,A V ) C\U,A V ), 
Ir f(z) := x J* z f(y)dy on C°(U,A V ) exists so that 

(1) ri?f(z) = f(z) 

for each z G U; this anti- derivative is H-linear (or H-left-linear when v = 2): 

(2) X r [af{z) + bg(z)} = al r f(z) + bX r g{z) 

for any real constants a,b G R (or a,b G H for v = 2) and continuous 
functions f, g G C°(U, A v ). If there is a second anti- derivative Ty^fiz), 
then Xyf(z) — Xy^f{z) belongs to the kernel ker(T) of the operator T. 

Proof. Using the multiplication on the marked doubling generator i^u 
from the right we have 

k-l k-1 

(3) S2ij(9g(z)/dz j )]i 2 «ip3(z) = ^{dg{z) j 'dz^ij+^jiz), 
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where iji^v =: ij+2 v for each < j < 2" — 1, 2 < v. On the other hand, 
E^i j (dg(z)/dz j )^ j (z) = [E k J =o(dg(z)/dz j yz*^(z)}*, since fy(z) is real 
for each j and z. Therefore, it is sufficient to consider the first-order partial 
differential operator of the form: 

(4) Tg(z) = Y? =1 {dg/dz j )i$1>j(z) 
on the R-linear space C l (U, A v ) of all continuously differentiable functions 
g : U — >■ A v by real variables zq, z^a-i, where < n < 2 V — 1. The space 
of super-differentiable functions is everywhere dense in C°(U, A v ) and the 
line integral has the continuous extension on C°(U, A v ) along any continuous 
rectifiable path in U. Therefore, we take the space of super-differentiable 
functions and then take the continuous extension of Zy on C°(U,A V ) such 
that 

lim; l T f l = I T linn f =X r f 
for a sequence / of super-differentiable functions uniformly converging to 
/ on compact sub-domains V in U, where Tyf 1 is described below. Each 
function ipj{ z ) is continuous and each function Vj(z) is continuously differ- 
entiable on U (see also below), consequently, the integral J ~ f l (y)d,Uj(y) is 
continuously differentiable by z = 7(1) (i.e. by each real variable z^) and 
their sequence by I uniformly converges on each compact sub-domain V in 
U. Therefore, from TXyf 1 = f for each natural number leNwe get 

Tl r f = TX T lim, /' = lim, Tl r f = lim, /' = /, 
since the sequence {Iyf l (z)\ v ■ 1} is fundamental in C l {V,Av) for each 
compact sub-domain V in U and Xyf(z) G C 1 (?7, A v ). 

Consider the left algorithm of a calculation of the line integral over the 
Cayley-Dickson algebra A v (see §3 and references therein). We shall seek an 
anti-derivative in the form: 

(5) t f f{y)dy := n- 1 j Z f(y)dvM)h 

and use the homotopy theorem in the domain U satisfying conditions 1(D1, D2) 
so that 7 is a continuous rectifiable path joining points qZ = 7(0) and 
z = 7(1) (see [22l E2J EO])- Moreover, a branch of the anti-derivative op- 
erator Xyf(z) can be chosen such that it can be expressed with the help of 
a non-commutative line integral. 
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In view of Theorem 2.11 |27] and §3 we get 

(6) (d(J ( ; z f(y)du ] (y))/dz k ) = Cf{z)\dv 3 {z)/dz k ]) 

(see also the chain rule over the Cayley-Dickson algebra in [22| 123} 120]). 

Next we need some identities in the Cayley-Dickson algebra. Each Cayley- 
Dickson number has the decomposition: z = z i + ... + -^"-l^V-ij where 
z , Z2«-i G R, 2 G A v . To establish the identity 

(7) (ay)z* + {az)y* = a2Re(yz*) 

for any a,y,z G A v it is sufficient to prove it for any three basic generators of 
the Cayley-Dickson algebra A v , since the real field R is its center, while the 
multiplication in A v is distributive (a+y)z = az + yz and ((aa)((3y))(-yz*) = 
(a(3'y)((ay)z*) for all a, /3,7 G R and a,y,z G A v . If a = i , then (7) is 
evident, since yz* + zy* = yz* + (yz*)* = 2Re(yz*). If either y — i , then 
(ay)z* + (az)y* = az* + az = a2Re(z) = a2Re(yz*). Analogously for z = %q. 
For three purely imaginary generators i p ,i s ,ik consider the minimal Cayley- 
Dickson algebra $ = alg-R,{ip^i s -,ik) over the real field generated by them. 
If it is associative, then it is isomorphic with either the complex field C 
or the quaternion skew field H, so that (ay)z* + (az)y* = a(yz* + zy*) = 
a2Re(yz*). If the algebra $ is isomorphic with the octonion algebra, then 
we use Formulas 2 (Ml, M2) for either a, y G H and z = 1 or a, z G H and 
y = 1. This gives (7) in all cases, since the algebra algn{i pi i s ) with two basic 
generators i p and i s is always associative. Particularly, if y = i s ^ z = if., 
then the result is zero. 

Using (7) we get more generally, that 

(8) ((ay)z*)b* + ((az)y*)b* = (a2Re(yz*))b* = (ab*)2Re(yz*), 
consequently, 

(9) ((ay)z*)b* + ((az)y*)b* + {{by)z*)a* + ({bz)y*)a* = ARe(ab*)Re(yz*) 
for any Cayley-Dickson numbers a, b,y,z G A v . 

We shall take unknown functions Vj(z) G A v as solutions of the system 
of linear partial differential equations by real variables z k . 

(10) duj(z)/dzj = l/ipj(z) for all 1 < j < n and z G U; 

(11) ipk(z)di>j(z) I 'dzk = ipj(z)di>k(z) I 'dzj for all 1 < j < k < n and z G U . 
Each function Uj(z) can be written as fj(z) = YhLq 1 v j,i{ z )ii with real- valued 
components Uji(z). Practically, it is sufficient to consider non-zero Uj t i(z) for 
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I = 1, n. Thus using the generators i , the system can be written 
in the real form. This system has a non-trivial C 1 solution Vj{z) for each j 
(see §12.2 [29], particularly, in the class of super-differentiable functions for 
super-differentiable ipj(z) see also [26, 20J). In System (10,11) functions ipj 
are real and coordinates are real, consequently, a solution {vj(z) : j} may 
be chosen real- valued. 

From Identities 3(2, 3) and (6, 9 — 11) we infer that 

(12) £ [(d( f f(y)du j (y))/dz k )i j ]il^) = 
{[(/(z).(%(^)/^ fe ))i i ]i^ fe (^)+[(/>)-(^^)/%))4]^(^)} = and 

l<j<k<n 

(13) f f(y)dv j (y))/dz j )i j ]i^ j (z) = nf(z), 

since Z)?=i — n an d n is some fixed natural number for the domain U, 
f(z).x = f(z)x for each real number x, (zij)i* = z for each z G A v , where / 
is the operator corresponding to dn(z)/dz, when / is in the ^-representation 
fi (see the notation in §3). Using Formulas (4, 5, 12, 13) we get Formula (1). 

From the identity / a\idz = a J fidz for a suitable branch of the line 
integral given by the left algorithm and for each non-trivial phrase \i and 
constants a, b G R for v > 3 or a, b G H for v = 2 (see the rules in [221 [231 
EH El]) we get Formula (2). 

Since T(Ty f(z)—Ir, 2 f(z)) = 0, the difference {Xyf{z)—Xt^f{z)) belongs 
to the kernel ker{T) = T^O), where T : C l (U,A v ) ->• C°(U,A V ). 

4.1. Example. If ipj depends only on Zj for each j, there exists a C l dif- 
ferentiable change of variables £ = ((z) so that dg(Q/dQ = (dg(z) / dzj)if)j(z) 
for each differentiable function g : U — )■ ^ by real variables £0, ^2«-i 011 
U, where 

(1) (dz k /d(j) = 5j !k ip k (z) 

for all j and k, djj = 1, while 5j± = for each j 7^ fc. We take new functions 
jg satisfying the equation: 

(2) jg(ijZ) = g(z) for each z G U and all j. We also put 

(3) Vj (z)=t*z. 

The multiplication of generators implies that i*j{ijz) = z for all j = 0, 2 V — 1 
and z G *4.„. Therefore, from Equations (1,2) we deduce that 
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(4) (dg(z)/dz).i j = (d j g(r] j )/dr] j ).\(dr lj /dz).i j \ = (d^fa) / 'c%).l = 
since (d?]j/dz).ij = = 1 for each j. Then we take the integral 

rz n rz 

(5) r g(y)dy.= n- l J2 j9(Vj(y))ijdVj(y), 

since Iozj9(vj(y))ijdVj(y) = (I z zj9(vj(y))dvj(y))ij- 

Mention that generally T(f(z)b) may be not equal to (Tf(z))b for a 
constant b G A v \ R and a function / G C 1 (U,A V ) with t> > 2, since the 
Cayley-Dickson algebra is non-commutative. 

This theorem can be generalized in the following manner encompassing 
wider class of partial differential operators of the first order over Cayley- 
Dickson algebras. 

5. Theorem. Suppose that the first order partial differential operator T 
is given by the formula 

(1) T/ = Y.U{df/d Zj )ct>*{z), 

where (f)j(z) ^ {0} for each z G U and (f)j(z) G C°(U, A v ) for each j = 0, ...,n 
such that Re((f>j(z)(f)l(z)) = for each z G U and each < j ^ k < n, where 
a domain U satisfies Conditions 1(D1,D2), z is a marked point in U, 
\<n<2 v ,2<v. Suppose also that the system {(p (z), (fi n (z)} is f or n = 
2 V — 1, or can be completed by Cayley-Dickson numbers <p n+ i(z), (p2 v -i(z), 
such that (a) algn{(j)j(z), (f>k{z), 4>i{z)} is alternative for all < j,k,l < 
2 V — 1 and (j3) algn{<f) (z), 4>2 v -i( z )} — A v for each z G U. Then a line 
integral X T : C°(U,A V ) C\U,A V ), l r f(z) := r f* z f(y)dy on C°(U,A V ) 
exists so that 

(2) Tlyf(z) = f(z) 

for each z G U; this anti- derivative is H-linear (orH-left-linear when v = 2). 
If there is a second anti- derivative X^^f{z), thenX~[f(z) — Ix,2f(z) belongs 
to the kernel ker(T) of the operator T. 

Proof. We shall demonstrate that a branch of the anti-derivative oper- 
ator X~cf(z) can be chosen such that it can be expressed with the help of 
a non-commutative line integral from §3. Using the technique of §4 we can 
consider the case of purely imaginary (f)j(z) for all z G U and j = 0, n. We 
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seek an anti-derivative operator in the form: 

(3) T / f(z)dz = (n + l)- 1 £ f 

For finding unknown functions q and i/j, j = 0, n we impose the following 
conditions: 

(4) (q(z) .[di/j(z) I 'dzj])(j)*(z) = f(z) for each j = 0, ...,n and 

(5) (g(z).[^<^/^ fc ])^(^) + (g(^).[^ fc (^/^])0*(^) = for all < j < 
k < n. 

As in §4 it is sufficient to consider the case of a locally analytic (super- 
differentiable) function / using the limit transition. The function / is given 
on U and it defines the operator / on U, i.e. its phrase jx is prescribed by the 
left algorithm for a given phrase fi of / (see [221 [231 [251 24J). The operator 
q means that a function g and a phrase k of g exist such that 
= dg(z)/dz, q(z). 1 = g(,z) for each z E U. 

In accordance with the conditions of this theorem the algebra alg-n((pj(z), 4>k( z )) 
is alternative for all < j < k < n and z £ U. Therefore, due to Condition 
(j3) Equations (4, 5) take the form: 

(6) (dg(z)/dz).[duj(z)/dzj] = f(z)(l/(j)*(z)) for each j = 0,...,n and 

(7) ((dg(z)/dz).[diy j (z)/dz k ])(f>t(z) + ((dg(z)/dz).[diy k (z)/dz j ])(j>*(z) = 
for all < j < k < n. 

Solutions of this system exist (see [261 120]). To be more concrete we 
impose additional relations: 

(8) dvj{z)/dzj = 4>j(z) for all j = 0, ....,n and z G U, 
consequently, the system of partial differential equations (6) becomes: 

(9) (dg(z)/dz).c/ )j (z) = /(z)(l/^(*)) for each j = 0, ...,n, 

since algn{<j)j(z), 4>k{ z ), 0z(^)} is alternative for all < j, k, I < 2 V — 1 and 
algn{(p (z), (p2v-i(z)} = A v for each z G U so that each Cayley-Dickson 
number £ G A v has the decomposition £ = £ o 0o("2) + ••• + ^-i^-i^) with 
real coefficients £ ; •••j£2«-i G R. 

Solving the latter system (9) one gets the function g(z) on U. Substituting 
the known function g in System (6,7) one gets a C 1 solution v§{z\...,v n {z) 
on U; or a super-differentiable solution, when (f>j(z) for each j and /(z) are 
super-differentiable on U. Mention that the function g depends R-linearly 
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on /, since the system of equations which was considered above is linear by 
/ and g. Thus the operator q depends R-linearly on /. 
Using Formulas (4, 5) and 4(6, 9) we deduce that 

(10) E [<9(f "q(y)dv s (y))/dz k ]<j>* k (z) = 

j^k>0 J ° z 

E {[q{z).{dv 3 {z)/dz k )W k {z) + [^.(du^/dzj)}^)} = and 

0<i<fc<n 

(11) £P( / q{y)dv j {y))/dz j W j {z) = ^[^.(du^/dz^z) = (n+l)f(z), 

3=0 J ° z j=0 

since Re{(f)j{z)(f)* k {z)) = for each z G U and each < j ^ k < n. 

The rest of the proof is analogous to that of Theorem 4. 

6. Corollary. Let suppositions of Theorem 5 be satisfied so that (f)j(z) = 
ou(z] ij)ipj(z) for each z G U , where u is an H-linear automorphism u> : A v — >■ 
A v mapping the standard base of generators {ij} into a base of generators 
{oj(z; ij) : j = 0, 2 V — 1}, \u>(z; ij)\ = 1, where ipj(z) satisfies conditions of 
theorem 4 for each j — 0, n. Then the first order differential operator 5(1) 
has an anti- derivative X T on C°(U,A V ). Two anti- derivatives of Theorems 4 
and 5 under these suppositions are related with the help of the automorphism 
uj . 

Proof. This follows immediately from Theorem 5. It remains to find 
a relation between two anti-derivatives for two different partial differential 
operators: 

(1) ^wf = T,]=o(9f/dzj)(f)*(z) 
and T given by equation 2(13). 

For each Cayley-Dickson number z = z i + ... + ^-i^-i G A v its image 
is io(y;z) = z N + z 1 N 1 + ... + z 2 v- 1 N 2 v- 1 , consequently oo(y; z*) = [u(y;z)\*, 
where Zj G R, Nj = Nj(y) := u(y;ij) for each j. Particularly N = i , 
since i§i 3 - = ij and u(y;ij) = cu(y;ioij) = u){y;io)uj{y;ij) for each j and y. 
Therefore, u(y;x) = x for each real number x G R, since u(y; 1) = 1 and 
the mapping u(y; *) is R-linear by the second argument, 1 = i Q . Therefore, 
applying the automorphism u we deduce that 

(2) T u f(z)= U (z;T 8 (z)), 

where co(z; s(z)) = f(z) for each z G U, that is s(z) = ui^iz; f(z)), co^iz; *) 
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denotes the inverse automorphism by the second argument for z G U. Let us 
take the function f(z) = x w Jo Z 9(y)dy, where g(z) is a continuous function. 
Then T 0J f(z) = g(z) for each z £ U and from (2) and 5(1, 2) one gets 

(3) u^iz-giz)) = Tw 2 _1 ( 2 i t w /„* g{y)dy) = T r J^ z U2 1 {y,g{y))dy, conse- 
quently, applying T / and u>(z; *) one also gets 

(4) r w J* z g{y)dy = u>(z; T j* z ^\y\ g{y))dy 

for each continuous function g on U . 

6.1. Remark. If in Theorem 5 drop Conditions (a,/3), then partial 
differential equations 5(4, 5) will be hard to resolve. 

To specify the anti-derivative operator Xy in Theorems 4 and 5 more 
concretely it is possible to choose a family of rectifiable continuous paths 
(or C 1 paths) {Y : z G U} such that 7 2 (0) = z and 7 2 (1) = z and 
lim 2 _M,sup re[0il] | 7 2 (r) - Y{t)\ = 0. 

Another more rigorous procedure is in providing a foliation of a domain 
U by locally rectifiable paths {7" : a G A}, where A is a set. We take 
for definiteness a canonical closed domain U in A v satisfying Conditions 
1(D1,D2). 

A path 7 :< a,b >— > U is called locally rectifiable, if it is rectifiable 
on each compact segment [c, e] C< a,b >, where < a,b >= [a,b] := {t G 
R : a < t < b} or < a, b >= [a, b) := {t G R : a < t < b} or < a, 6 >= 
(a, b] := {t G R : a < t < b} or < a,b >= (a, 6) := {t G R : a <t < b}. 

A domain [/ is called foliated by rectifiable paths {7" : a G A} if 7 :< 
a a , b a >— > U for each a and it satisfies the following three conditions: 

(Fl) UaeA7a(< a>a,b a >) = U and 

(F2) 7 a (< a Q , 6 a >) D 7^(< a^, 6^ >) = for each a ^ (3 e A. 
Moreover, if the boundary dU = cl(U) \ Int(U) of the domain U is non-void 
then 

(F3) du = (\J aeAl YM) U (U^a 2 7 /3 (&/3)), 
where A 1 = {a e A :< a a) bp >= [a a ,bp >}, A 2 = {a G A :< a a ,&g >=< 
0^,6/3]}. For the canonical closed subset U we have cl(U) = U = cl(Int(U)), 
where cl(U) denotes the closure of U in A v and Int(U) denotes the interior 
of U in A v . For convenience one can choose C 1 foliation, i.e. each 7° is of 
class C 1 . When U is with non-void boundary we choose a foliation family 
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such that UaeA7( a a) = dUi, where a set dU\ is open in the boundary dU 
and so that w\gui would be a sufficient initial condition to characterize a 
unique branch of an anti-derivative w =Iyf. 

When dU ^ a marked point qz can be chosen on the boundary dU and 
each point on the boundary can be joined by a rectifiable path in U with z. 
This foliation is justified by the formula: 

S 1 f(z)du(z) = J 7l f(z)dv(z) + J, 2 f(z)du(z) 
for each continuous function / on U and each function v of bounded varia- 
tion on U, for any rectifiable paths 7 1 : [01,61] — > U and 7 2 : [02,62] — > U 
so that a = a\ < b± = a<i < b 2 = b while 7 : [a, b] — > U is given piece- 
wise as 7(i) = 7 1 (£) for each t G [ai,6i] and ^(t) = 7 2 (t) for each t G 
[02,62]- Thus instead of f* z f{z)dv(z), i.e. / f(z)du(z) with 7(a) = z 
and 7(6) = z, we take J a < f(z)dv(z) for any [c, e] C< a a ,b a >. If 

/ I [c,ej 

lim c ^ aa ^ ba Xy«| f(z)dv(z) converges we denote it by J 7 « f(z)dv(z) and 

take instead of the family {j* 7 a| [ce] f{z)dv{z) : [c, e] C< a a , 6^ >}. Therefore, 

a branch of the anti-derivation operator prescribed by the family {(/ 7 a J2j q.{y)dvj{y)) '■ ol G 

A} or {(/ 7 a| [c e] Ej q{y)dvj{y)) '■ a G A; [c, e] C< a a , 6/3 >} is defined up to a 

function defined on the boundary dU when it is non-void or by convergence 

to a definite limit at infinity along paths, when U is unbounded in certain 

directions Hi] in the Cayley-Dickson algebra A v , rj G A v . 

Clearly, boundary conditions are necessary for specifying a concrete solu- 
tion or a branch of an anti-derivative, since in the definition of the line inte- 
gral / f{z)dv{z) the operator / is restricted to the condition f(z).l = f(z) 
for each z G U so it is defined up to a function of 2 V — 1 independent real 
variables (see also §3). In accordance with the formulas of §§4 and 5 the 
anti-derivation operators are defined up to functions of 2 V — 1 real variables 
after a suitable change of variables. For example, Y,]=o(dg(z) /dzj)i* = 
for g(z) = nzo + z\i\ + ... + z n i* w or jyj=o(dq(z) / 'dzj)i* = on the plane 
Zq — Z\ — ... — z n = for q(z) = z\ + z\i\ + ... + z\i* n . These functions can be 
written in the ^-representation due to Formulas 1(1 — 3). 

For concrete domains some concrete boundary conditions can be chosen 
(see also below). Mention, that a minimal necessary correct boundary con- 
ditions may be not on the entire boundary, but on its part. Otherwise, they 
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may be on some hyper-surface S in U of real dimension 2 V — 1 depending on 
the domain, for example, for an infinite cylinder C in both directions along 
its axis with S being the intersection of C with a hyper-plane perpendicular 
to its axis. 

Mention that the homotopy theorem for domains satisfying Conditions 
1(-D1,.D2) is accomplished for super-differentiable functions on U (see [23], 
[22]), but for a continuous function / on U it may certainly be not true. 
This is caused by several reasons. If a family of locally analytic functions 
f n converges to / uniformly on a compact sub-domain V in U a radius r™ 
of local convergence of a power series of f n in a neighborhood of a point 
x G V may tend to zero with n tending to the infinity. Phrases /i n in the 
^-representation corresponding to f n may be inconsistent on the intersection 
V x n V y of open neighborhoods V x and V y of different points x, y G V , when 
Vx H V y 7^ . Functions f n or their phrases p n may be with branching points 
in the domain U. That is functions f n accomplishing the approximation of 
/ may have several branches on U and a slit of U by a 2 V — 1 dimensional 
sub-manifold S n over R may be necessary to specify branches of f n . But the 
family S n with different n may be inconsistent and S n may depend of n. 

For super-differentiable functions f n operator valued functions f n are also 
super-differentiable. If / is only continuous non super-differentiable function 
on the domain U, then the operator valued function / is defined only in the 
sense of distributions [/,7; v) = L f{z)dv{z) for any rectifiable path 7 in U 
and each function v of bounded variation on U. Moreover, the homotopy 
theorem may be non true for generalized functions (see below). 

7. Particular case. We consider a phrase v which can be presented as 

(-P3) v = p(/i) with a right ^-linear (super)-differentiable phrase /1 and 
a projection operator p being an R-linear combination of the projection op- 
erators TTj. Particularly, p may be the identity operator or one of the ttj. 

For any z-differentiable phrase ip and constants a, b G A v we have / a(ip(z)b)dz = 
a ((/ 7 ip(z)dz)b) and J (a^(z))bdz = (a(J i/)(z)dz))b. Then in view of the ho- 
motopy theorem [231 122] Equation 3(2) implies for any such v = p(p) that 

(1) / T(y(z))dz = p([ [dp(z)/dz}.{J2i(dz/dz).t 3 }(i*^(z)dz)}) = 
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p( / [dfi(z) / dz].^2{ij(i*ijjj(z))dz}) — p( [dp,{z) / dz].[a{z)dz\) 




= p(n(z a (P))) - p(p(z a {a))) = u(za(P)) - u(z a (a)), 





where t = J2j ipj ■ 

For non right ^4„-linear z-differentiable phrase p Formulas (1 — 3) may 
already be not valid. Certainly common line integrals of z-differentiable 
phrases (functions) can be calculated by the general algorithms (see [221 [221 
EH [27]). A result of the line integration along a rectifiable path 7 in the 
domain U we denote as the composition of two functions 



where A and £ are two 2-differentiable functions on their domains, 7(0) = a, 
7(1) = p. Frequently one can use a Cayley-Dickson subalgebra Q isomorphic 
with either the quaternion skew field H or the octonion algebra O so that 
7(1) — 7(0) G Q and use the homotopy theorem. On the other hand, each 
rectifiable continuous path 7 in the domain U in the Cayley-Dickson algebra 
A v can be presented as a uniform limit of rectifiable continuous paths 7™ in 
U composed of segments parallel to axes Ri^, k = 0, 2 V — 1. Therefore, 



for any continuous function on U (see [231 122] ) . The functions A and £ depend 
on tp so in more details we denote them by A = A^, and £ = 
Thus the general integral of Equation 2(12) is: 



where Im(z) := z — Re(z), Re(z) := (z + z*)/2. The term 0y(im(£)) takes 
into account the non-commutativity for 2 < v and non-associativity for 3 < v 



(5) ]T / [{dp{z)/dz).i j ]\i*ij j {z)}dzo = I {dv{£)/dt).dt 




A(£(/?))-A(£(a)), 




(6) A*(&(a;)) 
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of the Cayley-Dickson algebra A v , since its center is the real field R = Z(A V ) 
for any v > 2. There is the bijective correspondence between A^,(^) and / 
which will be specified below. 

8. Transformation of the first order partial differential operator 
over the Cayley-Dickson algebras. 

To simplify the operator T and its particular variant a one can use a 
change of variables. We consider this operator in the form: 

(1) T/ = YgJidf/dzMz), 

with either rjj(z) = i*ipj(z) or rjj(z) = (p*(z) G A v for each j (see Theorems 
4 and 5 above). For it we seek the change of variables x = x(z) so that 

(2) T.f^\dx l /dz j )u j (z)=t l , 

where £; G A v is a constant for each I, for rjj not being identically zero, while 
u>j is chosen arbitrarily also z-differentiable so that the resulting matrix Q will 
not be degenerate, i.e. its rows are real-independent as vectors (see below). 
Certainly (dxi/dzj) G R are real partial derivatives, since x\ and Zj are real 
coordinates. We suppose that the functions r)j(z) are linearly independent 
over the real field for each z in the domain U. Using the standard basis of 
generators {ij : j — 0, 2 V — 1} of the Cayley-Dickson algebra A v and the 
decompositions Uj = Y^k u j,kik an d tj = J^k^j^k with real elements Uj^ and 
tj t k for all j and k we rewrite System (2) in the matrix form: 

(3) {dxi/dZj)i ) j =0t ... t 2v-iSl = T, 

where = (uJj : k)j,k=o,. ..,2 v -i, j,k=o,...,2 v -i- Suppose that the functions 

u)j(z) are arranged into the family {uj : j = 0, 2 V — 1} as above and are 
such that the matrix Q(z) is non-degenerate for all z in the domain U. For 
example, this is always the case, when |wj(-2)| > and Re[u)j(z)u>k(z)*] = 
for each j ^ k for each z G U. Here particularly ujj(z) = rjj(z) can also be 
taken for all j = 0, ...,2 V — 1 and z G U. Therefore, Equality (3) becomes 
equivalent to 

(4) (dx l /dz j ) ld=0 ,..., 2 ^ 1 = Til- 1 . 

We take the real matrix T of the same rank as the real matrix {u)j t k)j,k=o,...,2 v -i- 
Thus (4) is the linear system of partial differential equations of the first order 
over the real field. It can be solved by the standard methods |29j . 

We remind how each linear partial differential equation (3) or (4) can be 
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resolved. Write it in the form: 

(5) X ± (xi, ...,x n ,u)du/dx 1 +...+X n (x 1 , ...,x n ,u)du/dx n = R(x 1 , ...,x n ,u) 
with u and x±, ...,x n here instead of x\ and z , ...,Z2v-i in (3) seeking simul- 
taneously suitable R corresponding to ti^. A function u = u(xi, ...,x n ) de- 
fined and continuous with its partial derivatives du/dxi, ...,du/dx n in some 
domain V of variables xi,...,x n in R n making (5) the identity is called a 
solution of this linear equation. If R = identically, then the equation is 
called homogeneous. A solution u = const of the homogeneous equation 

(6) Xi(xi, x n , u)du/dx\ + ... + X n (xi, x n , u)du/dx n = 
is called trivial. Then one composes the equations: 

(7) dx 1 /X 1 (x) = dx 2 /X 2 {x) = ... = dx n /X n (x), 

where x = (xi,...,x n ). This system is called the system of ordinary differ- 
ential equations in the symmetric form corresponding to the homogeneous 
linear equation in partial derivatives. It is supposed that the coefficients 
Xi, ...,X n are defined and continuous together with their first order partial 
derivatives by xi, x n and that X ± , ...,X n are not simultaneously zero in a 
neighborhood of some point x°. Such point x° is called non singular. For 
example when the function X n is non-zero System (7) can be written as: 

(8) dxi/dx n = X 1 /X n ,...,dx n - 1 /dx n = X n _ x jX n . 

This system satisfies conditions of the theorem about an existence of integrals 
of the normal system. A system of n differential equations 

(9) dy k /dx = f k (x,yi, ...,y n ), k = 1, ...,n, 

is called normal of the n-th order. It is called linear if all functions 
depend linearly on yi,...,y n . Any family of functions yi,...,y n satisfying 
(9) in some interval (a, b) is called its solution. A function g(x, yi, y n ) 
different from a constant identically and differentiable in a domain D and 
such that its partial derivatives dg/dyi,...,d/dy n are not simultaneously zero 
in D is called an integral of System (9) in D if the complete differential 
dg = (dg/dx)dx + (dg / dy^dyi + ... + (dg / dy n )dy n becomes identically zero, 
when the differentials dy^ are substituted on their values from (9), that is 
(dg{x,y)/dx) + (dg/dy^f^x^y) + ... + (dg(x,y)/dy n )f n (x,y) = for each 
(x,y) G D, where y = (y±, ■■■,y n )- The equality g(x,y) = const is called the 
first integral of System (9). 
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It is supposed that each function fk(x, y) is continuous on D and satisfies 
the Lipschitz conditions by variables y\, ...,y n : 

(L) \fk(x,y) - fk(x,z)\ < C k \y - z\ 
for all (x, y) and (x, z) G D, where are constants. Then System (9) has 
exactly n independent integrals in some neighborhood D° of a marked point 
(x°,y°) in D, when the Jacobian d(g±, g n )/d(yi, ...,y n ) is not zero on D° 
(see Section 5.3.3 [29]). 

In accordance with Theorem 12.1,2 each integral of System (7) is a 
non-trivial solution of Equation (6) and vice versa each non-trivial solution 
of Equation (6) is an integral of (7). If g%(xi, x n ),...,g n -i(xi, ...,x n ) are 
independent integrals of (7), then the function 

(10) U = $(#!,...,£„_!), 

where $ is an arbitrary function continuously different iable by g±, ...,g n -i, is 
the solution of (6). Formula (10) is called the general solution of Equation 
(6). 

To the non-homogeneous Equation (5) the system 

(11) dx\/X\ = ... = dx n /X n = du/R 

is posed. System (11) gives n independent integrals g lf g n and the general 
solution 

(12) $(#i(xi, ...,x n ,u), ...,g n (xi, ...,x n ,u)) = 

of (5), where $ is any continuously differentiable function by gi, g n . If the 
latter equation is possible to resolve relative to u this gives the solution of 
(5) in the explicit form u = u(x\, ...,x n ) which generally depends on $ and 
gi, g n . Therefore, Formula (12) for different R and u and Xj corresponding 
to ti t k and Xi and Ujk respectively can be satisfied in (3) or (4), the variables 
Xj are used in (12) instead of Zj in (3, 4), where k — 0, 2 V — 1. 

Thus after the change of the variables the operator T takes the form: 

(13) Tf^Ef^idf/dx^ 

with constants tj G A v . Undoubtedly, also the operator T with j = 0, ...,n, 
2 V ~ 1 < n < 2 V — 1 instead of 2 V — 1 can also be reduced to the form 
Tf = J21=o(9f/dxj)tj, when the rank is rank(uj t k) = n + 1 in a basis of 
generators N , ...,N n , where Nq,...,Nsv_i is a generator basis of the Cayley- 
Dickson algebra A v . Particularly, if the rank is rank(ujj t k) = m < 2 V and T 
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has the unit upper left m x m block and zeros outside it, then tj = Nj for 
each j — 0, m — 1 can be chosen. 

One can mention that direct algorithms of Theorems 4 and 5 may be 
simpler for finding the anti-derivative operator Xx, than this preliminary 
transformation of the partial differential operator T to the standard form 
(13). 

9. Definitions. 

Let X and Y be two R linear normed spaces which are also left and right 
A r modules, where 1 < r. Let Y be complete relative to its norm. We put 
X® k := X ©r ... ©r X is the k times ordered tensor product over R of X. 
By L q ^ k (X® k , Y) we denote a family of all continuous k times R poly-linear 
and A r additive operators from X® k into Y. Then L q>k (X® k ,Y) is also a 
normed R linear and left and right A r module complete relative to its norm. 
In particular, L qj i(X,Y) is denoted also by L q (X,Y). 

We present X as the direct sum X = X i © ... © A^-i^-i, where 
X ,...,X 2 r-i are pairwise isomorphic real normed spaces. If A E L q (X,Y) 
and A(xb) = (Ax)b or A(bx) = b(Ax) for each x E X and 6 e «4 r , then an 
operator A we call right or left ^4 r -linear respectively. 

An R linear space of left (or right) k times A r poly-linear operators is 
denoted by L l<k (X® k ,Y) (or L r , k (X® k ,Y) respectively). 

As usually a support of a function g : S — > A r on a topological space S is 
by the definition supp(g) = cl{t e 5 : ^(i) 7^ 0}, where the closure is taken 
in S. 

We consider a space of test function T> := X>(R n , F) consisting of all 
infinite differentiable functions / : R n — > Y on R n with compact supports. 
A sequence of functions /„ G V tends to zero, if all /„ are zero outside 
some compact subset K in the Euclidean space R n , while on it for each 
k — 0, 1, 2, ... the sequence {f^ ■ n E N} converges to zero uniformly. Here 
as usually f^ k \t) denotes the k-th derivative of /, which is a k times R poly- 
linear symmetric operator from (R™)^ to Y, that is f^ k \t).(hi,...,h k ) = 

(t).(/i CT (i), /icr(fc)) E Y for each hi,...,h k E R™ and every transposition 
a : {1, ...,&;} — > {1, ...,&;}, a is an element of the symmetric group S k , t E 
R n . For convenience one puts = f. In particular, /^(^.(e^, ej k ) = 
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d k f{t)/dt jl ...dt jk for all 1 < ji,...,j k < n, where e j = (0, 0, 1, 0, 0) G 
R n with 1 on the j-th place. 

Such convergence in V defines closed subsets in this space V, their com- 
plements by the definition are open, that gives the topology on T>. The space 
V is R linear and right and left A r module. 

By a generalized function of class V := [T?(R n , Y)\' is called a continuous 
R-linear ^-additive function g : V — > A r . The set of all such functional 
is denoted by V. That is, g is continuous, if for each sequence /„ G V, 
converging to zero, a sequence of numbers g(f n ) —: [g, f n ) G A r converges to 
zero for n tending to the infinity. 

A generalized function g is zero on an open subset V in R n , if [g, f) =0 for 
each / G V equal to zero outside V. By a support of a generalized function 
g is called the family, denoted by supp(g), of all points t G R n such that 
in each neighborhood of each point t G supp(g) the functional g is different 
from zero. The addition of generalized functions g, h is given by the formula: 

(1) [g + h,f) :=[g,f) + [h,f). 

The multiplication g G T>' on an infinite differentiable function w is given 
by the equality: 

(2) [gw, f) = [g, wf) either for w : R" — > A r and each test function / G T> 
with a real image /(R n ) C R, where R is embedded into Y; or w : R n — > R 
and / : R n ->■ Y. 

A generalized function g' prescribed by the equation: 

(3) [g 1 , f) := —[g, /') is called a derivative g' of a generalized function g, 
where /' G X>(R n , L q (R n , Y)), g' G [D(R n , L q (R n , Y))]'. 

Another space B := B(R n , Y) of test functions consists of all infinite dif- 
ferentiable functions / : R n — > Y such that the limit lim| t |^ +00 \t\ m f^\t) = 
exists for each m — 0, 1, 2, j — 0, 1, 2, .... A sequence f n G B is called con- 
verging to zero, if the sequence \t\ m f^\t) converges to zero uniformly on 
R n \ .B(R n ,0, R) for each m,j = 0, 1,2, ... and each < R < +oo, where 
B(Z,z,R) := {y G Z : p(y,z) < R} denotes a ball with center at z of 
radius R in a metric space Z with a metric p. The family of all R-linear and 
^4 r -additive functionals on B is denoted by B' . 

In particular we can take X = A 1 ?, Y = Af with 1 < a, j3 G Z. Analo- 
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gously spaces V(U, Y), [V(U, Y)}', B(U, Y) and [B(U, Y)\ are denned for do- 
mains U in R™. For definiteness we write B(U, Y) = {f\ v : / G B(R n ,Y)} 
tmdV(U,Y) = {f\u: feV(R n ,Y)}. 

A function g : U — >■ A v is called locally integrable, if it is absolutely 
integrable on each bounded A measurable sub-domain V in U, i.e. 

j v \g(z) | A (dz) < oo, where A denotes the Lebesgue measure on U. 

A generalized function / is called regular if locally integrable functions 
j,kf l -,if 2 '■ U — >■ A v exist such that 

L/» = j u {J2 j ,k,ij,kf 1 {z) k uj{z) j f 2 (z)} q{3) X{dz), 
for each test function ui G B(U, Y) or u G V(U, Y) correspondingly, where 
uj = (iw, poo), ku(z) G A v for each z G U and all k, g(3) is a vector 
indicating on an order of the multiplication in the curled brackets and it 
may depend on the indices j, I — 1, a, k = 1, (3. 

We supply the space £>(R n , Y) with the countable family of semi-norms 

(4) p a:k (f) :=sup xeRn \(l + \x\) k d a f(x)\ 

inducing its topology, where k — 0, 1, 2, a = (cti, a n ), < a j G Z. On 
this space we take the space S'(R n , F)/ of all Y valued continuous generalized 
functions (functionals) of the form 

(5) / = M) + "- + /2"-i«2"-i and g = g i + ... + ^-i«2"-i, where and 
gj G i3'(R",F), with restrictions on £?(R",R) being real- or Ci = R © iR- 
valued generalized functions /o, fw-i, 9o, 92»-i respectively. Let = 
0o^o + ••• + (^-i^-i with O , 2 "-i G i3(R n , R), then 

(6) [/, 0) = Yl 2 kj2o[fj-,4'k)ikij- Let their convolution be defined in accor- 
dance with the formula: 

(7) [/ *g,<l>) = Ef fe =o([/j * 0*, 
for each G i3(R™, F). Particularly, 

(8) (/ * g){x) = f(x -y)* g(y) = f(y) * g(x - y) 

for all x, y G R™ due to (7), since the latter equality is satisfied for each pair 
fj and g k . 

10. The decomposition theorem of partial differential operators 
over the Cayley-Dickson algebras. 
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We consider a partial differential operator of order u: 
(1) Af(x) = £ a a (x)d a f(x), 

\a\<u 

where d a f = <9' a ' f(x) / 'dx^ ...dx^ n , x = x i + ...x n i n , Xj G R for each j, 
l<n = 2 r — 1, a = (ao, « n ), | cu | = a + ... + a n , < ctj G Z. By the 
definition this means that the principal symbol 

(2) A := J2 ^(x)d a 

\a\=u 

has a so that |a| = u and a Q (:r) G A r is not identically zero on a domain 
[/ in A r . As usually C k (U, A r ) denotes the space of k times continuously 
differentiate functions by all real variables xo, x n on U with values in A r , 
while the x-differentiability corresponds to the super-differentiability by the 
Cayley-Dickson variable x. 

Speaking about locally constant or locally differentiable coefficients we 
shall undermine that a domain U is the union of sub-domains W satisfying 
conditions 15(Dl,i - vii) and U j fl U k = dU j n <9f/ fe for each j ^ k. All 
coefficients a Q are either constant or differentiable of the same class on each 
Int(U j ) with the continuous extensions on UK More generally it is up to a 
C u or rr-differentiable diffeomorphism of U respectively. 

If an operator A is of the odd order u — 2s — 1, then an operator E of 
the even order u + 1 = 2s by variables (t, x) exists so that 

(3) Eg(t,x)\ t=0 = Ag(0,x) for any g G C u+1 ([c,d] x U,A r ), where t G 
[c, d] C R, c < < d, for example, Eg(t,x) = d(tAg(t,x))/dt. 

Therefore, it remains the case of the operator A of the even order u = 2s. 
Take z = zqIq + ... + ^"-i^-i G A v , Zj G R. Operators depending on a less 
set z^, z\ n of variables can be considered as restrictions of operators by all 
variables on spaces of functions constant by variables z s with s {/i, /„}. 

Theorem. Let A = A u be a partial differential operator of an even order 
u = 2s with locally constant or variable (locally) C s ' or x- differentiable on U 
coefficients a. a (x) G A r such that it has the form 

(4) Af = c Ut i{B Ut if) + ... + c u , k (B u>k f), where each 

(5) B UjP = B u , p fl + Q u -i,p 
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is a partial differential operator by variables x mul+ ... +mup _ 1+1 ,...,x mu il+ ... +mitiP 
and of the order u, m Ui0 = 0, c Uj k(x) G A r for each k, its principal part 

(6) B UiPfi = J2\a\=s &p,2a(x)d 2a 

is elliptic with real coefficients SL Pt2a (x) > 0, either < r < 3 and f G 
C U {U, A r ), or r > 4 and f G C U (U, R). T/jen £/jree partial differential opera- 
tors T s and T{ and Q of orders s and p with p < u — 1 with locally constant 
or variable (locally) C s ' or x- differ entiable correspondingly on U coefficients 
with values in A v exist and coefficients of the third operator Q may be gen- 
eralized functions, when coefficients of A are discontinuous locally constant 
or C s ' discontinuous on the entire U or when s' < s, r < v, such that 

(7) Af = T s {T\f) + Qf. 

Proof. Certainly we have ordQ u _ 1:P < u — 1, ord(A — A ) < u — 1. We 
choose the following operators: 

(8) T70r) = E E (d a f(x))[w* p ij p , a ] and 

P ~ 1 H<«> a q =0\/q<(m Ui i+...+m UiP -i+l) and q>(m Uil + ...+m UiP ) 

(9) rtf(x) = j: E (d a f(x))[w p r p , a i 

P ~ 1 H<«, ag=OVg<(m,„ i i+...+m UiP _i + l) and q>(m u> i + ...+m u>p ) 

where w 2 = c UjP for all p and ip Pj0l (x) = — a Pi2a (a;) for each p and x, w p G *4. r , 
i> P ,a{x) G ^4 r ^ and ip p>a (x) is purely imaginary for a Pj2a (^) > for all a and x, 
Re(w p Im(i(j Pia )) = for all p and a, Im(x) = (x — x*)/2, v > r. Here A r , v = 
A v /A r is the real quotient algebra. The algebra A r , v has the generators ij 2 r, 
j = 0, 2 v ~ r -l. A natural number v so that 2 v ~ r -l > E" - ( mp+q ~ l 



-ip=l ^q=0 \ q 

is sufficient, where (^j = m\/(q\(m — g)!) denotes the binomial coefficient, 
(m+q-i^j - g ^ e number of different solutions of the equation a\ + ... + a m = q 
in no n- negative integers a,-. We have either d a+/3 f G Ar for < r < 3 
or d a+l3 f G R for r > 4. Therefore, we can take ip Pj0l (x) G i^qR, where 
q = q(p, a) > 1, q(p\ a 1 ) ^ q(p, a) when (p, a) ^ (p 1 , a 1 ). 

Thus Decomposition (7) is valid due to the following. For b = d a+l3 f(z) 
and 1 = i 2 r p and wGi, one has the identities: 

(10) (b(w\))(w*\) = ((wb)\)(w*\) = -w(wb) = -w 2 b and 

(11) (((bl)w*)l)w = (((bw)\)l)w = —(bw)w = —bw 2 in the considered 
here cases, since Ar is alternative for r < 3 while R is the center of the 
Cayley-Dickson algebra (see Formulas 2(M1,M2)). 
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This decomposition of the operator A 2s is generally up to a partial dif- 
ferential operator of order not greater, than (2s — 1): 

(12) Qf(x) = J2U c u,pQu-i, P + 

E|«|< S ,|,|< S;7 < a ,e< /3 ,| 7+ e|>0[n?:0 1 ft) (g)K^ + ^-7(*)) 

\(d^ a (x))((d%(x)}, 

where operators X s and X* are already written in accordance with the general 
form 

(13) T*f(x) = E\ al < s (d«f(x))r] a (x); 

(14) TJ/(x) = E ms (d^f(x)) V }(x). 

The coefficients of Q may be generalized functions, since they are calcu- 
lated with the participation of partial derivatives of the coefficients of the 
operator T*, but the coefficients of the operators X s and X* may be discon- 
tinuous locally constant or C s ' discontinuous on the entire U or s' < s when 
for the initial operator A they are such. 

When A in (3) is with constant coefficients, then the coefficients w p and 
tp Pi a for T rn and X™ can also be chosen constant and Q — J2p=i c u , P Qu-i,p — 0. 

11. Corollary. Let suppositions of Theorem 10 be satisfied. Then a 
change of variables locally affine or variable C 1 or x-differentiable on U cor- 
respondingly on U exists so that the principal part A 20 of A 2 becomes with 
constant coefficients, when a Pj2Q , > for each p, a and x. 

12. Corollary. If two operators E = A 2s and A = A 2s _i are related by 
Equation 10(3), and A 2s is presented in accordance with Formulas 10(4,5), 
then three operators T s , T 5 " 1 and Q of orders s, s — 1 and 2s — 2 exist so 
that 

(1) A 2s _ 1 = T s T s - 1 + Q. 

Proof. It remains to verify that ord(Q) < 2s — 2 in the case of A 2s -i, 
where Q = {d(tA 2s ^i)/dt — T s T*}| t=0 . Indeed, the form \(E) corresponding 
to E is of degree 2s — 1 by a; and each addendum of degree 2s in it is of 
degree not less than 1 by t, consequently, the product of forms A(T S )A(TJ) 
corresponding to X s and X* is also of degree 2s — 1 by x and each addendum 
of degree 2s in it is of degree not less than 1 by t. But the principal parts of 
\(E) and A(X S )A(XJ) coincide identically by variables (t,x), hence ord({E — 
T s Tl}\ t=0 ) <2s — 2. Let a(t,x) and h(t,x) be coefficients from X* and X s . 
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Using the identities 

a(t,x)d t d 1 tg(x) = a(t ) x)d 1 g(x) and 

h(t,x)d /3 d t [a(t,x)d^g(x)] = h(t,x)d^[(d t a(t,x))d^g(x)\ 
for any functions g{x) G C 2s ~ l and a(t,x) G C s , 

ord[(h{t,x)dP), (a{t,x)d*)]\ t=0 < 2s - 2, 
where d t = d/dt, \p\ < s - 1, \j\ < s, [A, B] := AB - BA denotes the 
commutator of two operators, we reduce (Y S T\ + Q\)\t=o from Formula 10(7) 
to the form prescribes by equation (1). 

13. We consider operators of the form: 

(1) (T k + pi r )f(z) := {Eo<\a\<k(d a f(z) Va (z)} + f(z)P(z), 

with r) a (z) G A v , a = (a , a^-i), < aj G N for each j, \a\ — a + ... + 

«2-l, PI r f(z) ■= f(z)P, 

&*f(z) := d\ a \f{z)/dz^...dz«?_- l \ 2<r<v<oo, (3(z) G A v , z , z 2 r_, G 
R, z = z i + ... + Z2r-\iy-\. 

Proposition. The operator (T k + /3)*(T k + j3) is elliptic on the space 
C 2k (R? r , Av), where (Y k + f3)* denotes the adjoint operator (i.e. with adjoint 
coefficients). 

Proof. In view of Formulas (1) and 4(8) the form corresponding to the 
principal symbol of the operator (T fc + (3)*(T k + 0) is with real coefficients, 
of degree 2k and non-negative definite, consequently, the operator (T fc + 
(3)*(Y k + 0) is elliptic. 

14. Example. Let T* be the adjoint operator defined on differentiable 
A v valued functions / given by the formula: 

(1) (T + P)*f = [±(df\z)/dz 3 )<P 3 (z)} + f(z)f3(z)*. 

j=0 

Thus we can consider the operator 

(2) R fi := (T + 0)(T + fl\ 

From Proposition (13) we have that the operator is elliptic as classified by 
its principal symbol with real coefficients. Put S = S . In the x coordinates 
from §8 it has the simpler form: 

n 

(3) (T + /?)(T + /?)7 = E(W&3)l*;r 

j=0 

n n 

+2 (d 2 f/dx j dx k )Re(t j tl)+2j2(df/dx j )Re(tp)HfW\ 2 +^^ 

0<j<k<n j=0 j=0 
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because the coefficients tj are already constant. After a change of variables 
reducing the corresponding quadratic form to the sum of squares J2j e j s "j we 
get the formula: 

(4) TT*f = ET^f/ds])^, 

where Sj G R, €j = 1 for 1 < j < p and €j = —1 for each p < j < m, m < 2 V , 
1 < p < m depending on the signature (p, m — p). 
Generally (see Formula 5(1)) we have 

(5) A = (T + /3)(T 1 + p)f{z) = B f(z) + Qf(z), where 

(6) B f(z) = E^fW/d^WzTWkiz) + [f(z)P 1 (z)]P(z) and 

(7) Qf(z) = E^[(9/(z)/^)(^)(^)V^)]^(^)+E J [(9/(^)/^)0)(^)1/3(^) + 

ZMizXdpw/dzkMiz), 

(8) (Ti + (P)f(z) = [Y, j {df{z)/dz j )<t>){zr\ + f(z)0\z). 

The latter equations show that coefficients of the operator Q may be 
generalized functions, when <f)j(z) for some j or (3 l {z) are locally C° or C° 
or locally C 1 functions, while 4>k(z) for each k and (5(z) are locally C° or C° 
functions on U . We consider this in more details in the next section. 

15. Partial differential operators with generalized coefficients. 

Let an operator Q be given by Formula 14(7) on a domain U. Initially it 
is considered as a domain in the Cayley-Dickson algebra A v . But in the case 
when Q and / depend on smaller number of real coordinates z , ...,z n -i we 
can take the real shadow of U and its sub-domain V of variables (zq, z n -i), 
where Zk are marked for example being zero for all n < k < 2 V — 1. Thus we 
take a domain V which is a canonical closed subset in the Euclidean space 
R n 2 «-i < n < 2 V - 1, v > 2. 

A canonical closed subset P of the Euclidean space X = R™ is called a 
quadrant if it can be given by the condition P := {x G X : qj(x) > 0}, where 
(qj : j G Ap) are linearly independent elements of the topologically adjoint 
space X*. Here Ap C N (with card(Ap) = k < n) and k is called the index 
of P. If x G P and exactly j of the g^'s satisfy qi(x) = then x is called a 
corner of index j. 

That is P is affine diffeomorphic with P n = n™=i[a.j, fy], where — oo < 
3j < < oo, [aj,fej] := {x G R : 3j < x < bj} denotes the segment in R. 
This means that there exists a vector p G R n and a linear invertible mapping 
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C on R n so that C(P) - p = P n . We put P' 1 := (t u ...,tj, ...,t n : tj = a,), 
P' 2 := (ti, ...,tj, ...,t n : tj = bj). Consider t = (tx, ...,t n ) G P n . 

This permits to define a manifold M with corners. It is a metric sepa- 
rable space modelled on X = R n and is supposed to be of class C s , 1 < s. 
Charts on M are denoted (Ui,ui,Pi), that is, ui : U\ —> Ui(U{) C Pi is 
a C s -diffeomorphism for each /, TJ\ is open in M, ui o Uj^ 1 is of C s class of 
smoothness from the domain Uj(Uif]Uj) ^ onto ui(UiDUj), that is, Ujou^ 1 
and ui o uj 1 are bijective, Uj Uj = M. 

A point x G M is called a corner of index j if there exists a chart (U, u, P) 
of M with x G Z7 and is of index ind,M{x) = j in w(Z7) C P. A set of 
all corners of index j > 1 is called a border <9M of M, x is called an inner 
point of M if ind M (x) = 0, so <9M = \Jj>id j M, where 9»Af := {x G M : 
indu{x) = j} (see also [30J ) . We consider that 

(.Dl) V is a canonical closed subset in the Euclidean space R n , that is 
V = cl(Int(V)), where Int(V) denotes the interior of V and cl(V) denotes 
the closure of V. 

Particularly, the entire space R n may also be taken. 

Let a manifold W be satisfying the following conditions (i — v). 

(i) . The manifold W is continuous and piecewise C a , where C l denotes 
the family of I times continuously differentiable functions. This means by 
the definition that W as the manifold is of class C° fl C^ oc . That is W is of 
class C a on open subsets Woj in W and W \ (Uj Woj) has a codimension 
not less than one in W. 

(ii) . W = \J?= o Wj, where W = \J k W Q>k , Wj D W k = for each k ^ j, 
m = dimftW, dim R Wj = m — j, Wj+i C dWj. 

(Hi). Each Wj with j = 0, ...,m — 1 is an oriented C a -manifold, Wj is 
open in {Jk=j W k - An orientation of Wj+i is consistent with that of dWj for 
each j = 0, 1, m — 2. For j > the set is allowed to be void or non-void. 

(iv). A sequence W k of C a orientable manifolds embedded into R n , 
a > 1, exists such that W k uniformly converges to W on each compact 
subset in R n relative to the metric dist. 

For two subsets B and E in a metric space X with a metric p we put 

(1) dist(B, E) := max{sup bgB dist({b}, E), sup egE dist(B, {e})}, where 
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dist({b}, E) := mi e€E p(b, e), dist(B, {e}) := inf bgB p{b, e), b G B, e G E. 

Generally, dim^W = m < n. Let (e*(x), e^(x)) be a basis in the 
tangent space T x W k at x G W consistent with the orientation of W k , k G N. 

We suppose that the sequence of orientation frames (e* (xk), e^(xfc)) 
of PU fc at Xfc converges to (ei(x), ...,e m (x)) for each a; G Wo, where lim fc Xfc = 
x G Wo, while ei(x),...,e m (x) are linearly independent vectors in R n . 

{y). Let a sequence of Riemann volume elements A& on W k (see §XIII.2 
[55] ) induce a limit volume element A on IU, that is, X(BnW) = lim fc _ i . 00 (Sn 
IU fc ) for each compact canonical closed subset B in R n , consequently, \{W\ 
W$) = 0. We shall consider surface integrals of the second kind, i.e. by the 
oriented surface W (see (iv)), where each Wj, j = 0, m — 1 is oriented (see 
also §XIII.2.5 [SB]). 

Suppose that a boundary dU of U satisfies Conditions (i — v) and 

(vii) let the orientations of dll k and U k be consistent for each k G N (see 
Proposition 2 and Definition 3 [38J). 

Particularly, the Riemann volume element A& on dU k is consistent with 
the Lebesgue measure on U k induced from R n for each k. This induces the 
measure A on dU as in (y). This consideration certainly encompasses the 
case of a domain U with a C a boundary dU as well. 

Suppose that Ui,...,Ui are domains in R n satisfying conditions (Dl,i — 
vii) and such that Uj fl Uk = dUj fl dUk for each j ^ k, U = Uj=i Uj. 
Consider a function g : U — > Av such that each its restriction g\u is of 
class C s , but g on the entire domain U may be discontinuous or not C k , 
where < k < s, 1 < s. If x G dUj fl 9C/jt for some j ^ k, x G Int(U), 
such that x is of index m > 1 in Uj (and in Uk also). Then there exists 
canonical C a local coordinates (yi, ...,y n ) in a neighborhood 14 of x in {7 
such that S := V x n <9 m ^ = : y G V^; j/i = 0,...,y m = 0}. Using 
locally finite coverings of the locally compact topological space dUj fl dUk 
without loss of generality we suppose that C a functions Pi(z), P m (z) on 
R n exist with S = {z : z G R n ; P x (z) = 0, P m (z) = 0}. Therefore, on the 
surface S the delta-function 5(P\, P m ) exists, for m = 1 denoting them 
P — Pi and S(P) respectively (see §111.1 [6]). It is possible to choose yj = 
Pj for j = l,...,m. Using generalized functions with definite supports, for 
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example compact supports, it is possible without loss of generality consider 
that yx, y n G R are real variables. Let 9{Pj) be the characteristic function 
of the domain Vj := {z : Pj(z) > 0}, 9(Pj) := 1 for Pj > and 9(Pj) = 
for Pj < 0. Then the generalized function 9 (Pi, ■■■,P m ) '■= 9(Pi)...9(P m ) can 
be considered as the direct product of generalized functions 9(yi),...,9(y m ), 
l(y m+ i, y n ) = 1, since variables yi,...,y n are independent. Then in the 
class of generalized functions the following formulas are valid: 

(2) 89 (Pj) I 'dzk = 5(Pj)(8Pj/8z k ) for each k = 1, ...,n, consequently, 

(3) grad[9(P 1 ,...,P m )]=T l f =1 [9(Pi)..^P j -i)S(P j )(grad^ 

where , grad g(z) := (dg(z)/dz\, ...,8g(z)/8z n ) (see Formulas 111.1.3(1, 7, 7', 9) 
and 111.1.9(6) ©). 

Let for the domain U in the Euclidean space R n the set of internal surfaces 
clu[Int-R,n(U) fl [jj^k(8Uj fl 8Uk)} in U on which a function g : U — > A v or its 
derivatives may be discontinuous is presented as the disjoint union of surfaces 
Tj, where each surface T- 7 is the boundary of the sub-domain 

(4) V* := {P h i(z) > 0, P j>m .(z) > 0}, P = dV* = l)Zi d k V\ 

rrij e N for each j, clx(V) denotes the closure of a subset V in a topolog- 
ical space X, IntxiV) denotes the interior of V in X. By its construction 
{V^ : j} is the covering of U which is the refinement of the covering {11% : k}, 
i.e. for each "P- 7 a number k exists so that C Uf. and dV^ C dUk and 
Ui V j = UkU k = U. We put 

(5) hj(z(x)) = h(x)\ xeV i := 

'i" 1 //,:,!) .n g(z(x + y)) - liin, ; : .„ . g(z(x - y)) 

in accordance with the supposition made above that g can have only discon- 
tinuous of the first kind, i.e. the latter two limits exist on each T- 7 , where x 
and y are written in coordinates in V\ z = z(x) denotes the same point in 
the global coordinates z of the Euclidean space R n . We take new continuous 
function 

(6) g 1 (z)=g(z)-E j h j (z)9(Pj il (z),...,P j , mj (z)). 

Let the partial derivatives and the gradient of the function g 1 be calculated 
piecewise one each Uk- Since g 1 is the continuous function, it is the regular 
generalized function by the definition, consequently, its partial derivatives 
exist as the generalized functions. If g 1 \u j £ C x (Uj, A v ), then dg 1 (z) / dzk is 
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the continuous function on Uj , i.e. in this case dg 1 (z)xuj ( z ) I ® z k is the regular 
generalized function on Uj for each k, where Xg( z ) denotes the characteristic 
function of a subset G in A V) Xg( z ) = 1 for each z G G, while x( z ) — 
for z G A v \ G. Therefore, g 1 (z) = T,j 9 1 ( z )Xu j \\J k< .u k ( z ): where u o = 0, 

j, keN. 

On the other hand, the function g(z) is locally continuous on U, conse- 
quently, it defines the regular generalized function on the space V(U,A V ) of 
test functions as 

[g,u) := J u g(z)uj(z)X(dz), 
where A is the Lebesgue measure on U induced by the Lebesgue measure on 
the real shadow R 2 " of the Cayley-Dickson algebra A V) uj G V(U,A v ). Thus 
partial derivatives of g exist as generalized functions. 

In accordance with Formulas (2, 3, 6) we infer that the gradient of the 
function g(z) on the domain U is the following: 

(7) grad g(z) = grad g 1 (z) + J2j hj(z)grad 9(P j:1 , Pj, mj ). 

Thus Formulas (3, 7) permit calculations of coefficients of the partial dif- 
ferential operator Q given by Formula 14(7). 

16. Line generalized functions. 

Let U be a domain satisfying Conditions 1(D1,D2) and 15(Dl,i — vii). 
We embed the Euclidean space R n into the Cayley-Dickson algebra A v , 
2'"^ 1 < n < 2 V — 1, as the R affine sub-space putting R™ 3 x = (x± : ...,#„) H> 
Xii^ + ... + x n ij n + x° G A v , where jk ^ ji for each k ^ I, x° is a marked 
Cayley-Dickson number, for example, jk = k for each k, x° = 0. Moreover, 
each Zj can be written in the ^-representation in accordance with Formulas 
1(1-3). 

We denote by P = ~P(U) the family of all rectifiable paths 7 : [a 7 , 6 7 ] — > U 
supplied with the metric 

(1) p( 7 , u) := | 7 (a) - oj( aul ) \ + inf V^t)) - u(cf>(t)) 
where the infimum is taken by all diffeomorphisms : [a 7 , 6 7 ] — > [a u ,b u ] so 
that 0(a 7 ) = a w , a = a 7 < 6 7 = b (see §3). 

Let us introduce a continuous mapping g : B(U, A v ) x P(U) x V(U, A v ) — >■ 
Y such that its values are denoted by [g;u>, 7; v], where Y is a module over 
the Cayley-Dickson algebra A v , uj G B(U,A v ), 7 G P(i7), V(C/, ^) denotes 
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the family of all functions on U with values in the Cayley-Dickson alge- 
bra of bounded variation (see §3), v G V(U,A V ). For the identity mapping 
v{z) = id(z) = z values of this functional will be denoted shortly by [g; u, 7]. 
Suppose that this mapping g satisfies the following properties (Gl — G5): 

(Gl) [g; uj, 7; u] is bi- R homogeneous and A v additive by a test function 
uj and by a function of bounded variation v\ 

(G2) [g; u, 7; u] = [g; u, 7 1 ; v\ + [g; u, 7 2 ; u] for each 7, 7 1 and 7 2 G P(C/) 
such that 7(t) = 7 X (t) for all t G [a 7 i,6 7 i] and 7(t) = 7 2 (t) for any t G 
[a 7 2, 6 7 2] and a 7 i = a 7 and a 7 2 = 6 7 i and 6 7 = 6 7 2. 

(G3) If a rectifiable curve 7 does not intersect a support of a test function 
ou or a function of bounded variation v, j([a,b] fl (supp(u) fl supp(u)) = 0, 
then [g;o;,7; z/] = 0, where suppiu) := c/{^ G U : 7^ 0}. 

Further we put 

(G4) [flM^VflC-^r^.T] := (-i) l 1 5 ;9w w (z)/& m »...& 2 rr i l7 ] 

for each m = (m , ...,m 2 ^-i), is a non-negative integer < rrij G Z for 
each j, \m\ := m + ... +m 2 »_i. In the case of a super-differentiable function 
w and a generalized function <? we also put 

(G5) [(^(z)/^).^,...,^)^^] := (-1)% (d fc a;(^)/^ fc ).(/ii,...,/i fc ),7] 
for any natural number k G N and Cayley-Dickson numbers hi, ..,h^ G Ay. 

Then 5 is called the Y valued line generalized function on B(U, A v ) x 
P(C/) x V(U,Av). Analogously it can be defined on V(U,A V ) x P(U) x 
V(C/, A„) (see also §9). If Y — Av we call it simply the line generalized 
function, while for Y = L q (A^,A l v ) we call it the line generalized operator 
valued function, k, I > 1, omitting "on B(U, Ay) xP([/)x V(C/, A„)" or "line" 
for short, when it is specified. Their spaces we denote by L q (B(U,A v ) x 
~P(U) x V(U, Ay); Y). Thus if g is a generalized function, then Formula (G5) 
defines the operator valued generalized function d k g(z)/dz k with k G N and 
I = 1. 

If g is a continuous function on U (see §3), then the formula 

(G6) [0;w,7;i/] = J J g(y)uj(y)du(y) 
defines the generalized function. If f(z) is a continuous L q (A v ,A v ) valued 
function on U, then it defines the generalized operator valued function with 
Y = L q (A v ,A v ) such that 
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(G7) [f-,u, r ,v} = J 7 {f(z).uj(z)}dv(z). 
Particularly, for v = id we certainly have dv[z) = dz. 

We consider on L q (B(U, A v ) x P(U) x V(U,A V ); Y) the strong topology: 

(G8) lini; f l — f means that for each marked test function u G <B(C/, A) 
and rectifiable path 7 G P{U) and function of bounded variation v G V(£/, A) 
the limit relative to the norm in F exists 

limif/^w,^!/] = [/;w,7;i/]. 

17. Line integration of generalized functions. 

Let C^(V, A v ) denote the R linear space and right and left A v mod- 
ule of all functions 7 : V — > A v such that 7(2) and each its derivative 
d^g(z)/dz^ 11 ...dz™ n for 1 < < m is absolutely continuous on V (see 
§§3 and 16). This definition means that C m+1 (\/ A v ) C C^(V, A), where 
C m (V, Av) denotes the family of all m times continuously differentiable func- 
tions on a domain V either open or canonical closed in R™, which may be a 
a real shadow of U as well. 

17.1. Lemma. Let 7 G C^([a,o], A) n P(C/) and w G #([/, A) 
v G C® b (U, Av) for m = or u = id for m > 1, where < m G Z, £/ien a 
line generalized function [g; u>, 7|[ a ,x]i ^] continuous for m = or 0/ c/ass 
C m fry i/ie parameter x G [a, 6] /or m > 1. 

Proof. For absolutely continuous functions j(t) and z/ (i.e. when m — 0) 
the continuity by the parameter x follows from the definition of the line 
generalized function, since liniA^o p(l\[a,x], 7l[o,x+Ax]) = and 

lim A:E ^oP(^ 7l[a,z],^ 7l[a,x+A 2; ]) = 0. 

Consider now the case m > 1 and z/ = id. In view of properties 16(G1, G2) 
for any Ax 7^ so that x G (a, 6] := {£ G R : a < t < 0} and x + Ax G 
(a, 6) := {t G R : a < t < 6} the difference quotient satisfies the equalities: 

(1) {[g;u,i\[a,x+Ax]]-[g;uj,i\[a,x]\}/ Ax = id; u/Ax,-fO(j)\ [aiX] ]-[g; u/Ax,j\ [ajX] ], 

where : [a, x] — >■ [a, x + Ax] is a diffeomorphism of [a, x] onto [a, x + Ax] 
with 0(a) = a. Therefore, Acu := ou(z + Az) — u>(z) for z = j(t) and 
z + Az = 7(0(t)), t G [a, x] in the considered case. Using Conditions (Gl, G3) 
one can mention that if u = u l on an open neighborhood V of 7 in U, then 

(2) [ff;w,7] = [9;^\i], 

since w-w 1 = 0onV and 7 fl supp(u — u 1 ) = 0. 
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From Conditions 16(G1,G4) and Formula (2) we deduce that 

(3) limAa^oit^^l^+A*]] - [g;u,i\{a,x]}}/Ax 

= E^lo 1 ^; (du(z)/dzj), {dnfj(t)/dt)y\ M ], 
where z/ = djj(t)/dt for z = 7(i), t G [a, b], since each partial derivative 
of the test function u is again the test function. From the first part of the 
proof we get that [g; co, 7|[ ,2]] is of class C 1 by the parameter x, since the 
product (d , jj(t)/dt)'y(t) of absolutely continuous functions (d^j{t)/dt) and 
7(t) is absolutely continuous for each j. Repeating this proof by induction 
for k = 1, m one gets the statement of this lemma for 7 G C^([a, b],A v )C\ 
P(U). 

17.2. Lemma. If 7 is a rectifiable path, then a line generalized function 
[g; u, 7|[o,a;]] is of bounded variation by the parameter x G [a, b}. 

Proof. Let 7 G P(U) be a rectifiable path in U, 7 : [a, b] — > ?7. We can 
present 7 in the form 

(1) 7W = E?lo 1 7,(^, 

where each function jj(t) is real- valued. Therefore, jj(t) is continuous and 
of bounded variation for each j, since 7 is such. Thus the function u;(7(£)) 
is of bounded variation V^ 6 u;(7)) < 00, since oj is infinite different iable and 
7 ([a, 6]) is compact. 

On the other hand, each function / : [a, b] — > R of bounded variation can 
be written as the difference / = f 1 — f 2 of two monotone non- decreasing 
functions f 1 and f 2 of bounded variations: := V*f and / 2 (t) = — 

f(t) for each t G [a, 6] (see US])- This means that f k = g k + h k , where a 
function g fc is continuous monotone and of bounded variation, while h k is a 
monotone step function, where k — 1,2. When the function / is continuous 
one gets / = g 1 — g 2 . For a monotone non-decreasing function p one has 

VaP = P{t) - Pip). 

In view of Property 17(d) we infer that 

(2) [g;^n\ia,x]] =£S 1 [&sw,7|[ a , !B ]]*j, 

where the function [^-; a>, 7|[ a ,cc]] by x is real-valued for any uj G B(U, A v ) and 
7 G P(U) for all j = 0, ...,2 V - 1. 

The metric space P(t/) is complete, where U = cl{U). Indeed, let g n be a 
sequence of rectifiable paths in U fundamental relative to the metric p given 
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by Formula 16(1). Using diffeomorphism preserving orientations of segments 
we can consider without loss of generality that each path g n is defined on the 
unit segment [0, 1], a = 0, b = 1. It is lightly to mention that 

(3) \g(a) - f(a))\ + V a b (g - f) > sup teM \g(t) - f(t)\ 
for any two functions of bounded variation, /, g : [a, b] — > U . For each e > 
a natural number uq = no(e) exists so that p{g n ,g m ) < e/2 for all n,m> uq. 
That is (j) n : [0, 1] — > [0, 1] diffeomorphisms exist such that 

\g n (a) — g m (a))\ + V^(g n o <p n — g m o <p m ) < e for all n,m > n , since 
(f) m o is also the diffeomorphism preserving the orientation of the seg- 

ment. Using induction by e = l/l with I £ N one chooses a sequence of 
diffeomorphisms (f) n such that for each I £ N a natural number n Q = n (l) 
exists such that 

\g n (a) - g m (a))\ + V*(g n o f - f o m ) < l/l for all n,m > n (l), 
consequently, 

su PteM \9 n (<i> n (t)) - 9 m {4> m {t))\ < l/l for all n,m > n (l). 
Thus the sequence g n o <f> n is fundamental in C°([a, b], U). The latter metric 
space is complete relative to the metric 

d(f,g) := sup teM \f{t)-g(t)\, 
since from the completeness of the Cayley-Dickson algebra A v considered as 
the normed space over the real field the completeness of the closed subset U 
follows (see also Chapter 8 in [3]). Therefore, the sequence g n o n converges 
to a continuous function / : [a, b] — > U . On the other hand, linim^oo p(g n o 
( f) n J g m o m ) = p(g n o 0™, /) < l/l for each n > n (l), I £ N. The function 
g n ocf) n is of bounded variation, consequently, the function / is also of bounded 
variation. That is / £ ~P(U). Thus P(U) is complete. 

Take any sequence 7 n of C^ b ([a,b], A v ) paths in U converging to 7 rela- 
tive to the metric p on P(U) and the latter metric space is complete as it 
was demonstrated above. In view of Formula 17.1(3) and Property 16 (G3) 
the sequence [g; cj, 7 n |[ ,a;]] is fundamental in P(U). On the other hand, the 
generalized function g is continuous on B{U,A r ) x P(U), consequently, the 
sequence [g;u>, 7 n |[ a ,s]] converges in B(U, A r ) xP(U) to [g;uj, 7 n |[ a ,a;]] for each 
a < x < b, hence [g; u>, 7|[ ,a;]] = h m n[fl'; l n \[a,x}] in P(U). By the conditions 
of this lemma [g; u, 7|r a ,a;]] ^ P{U), since 7([a, b]) C U. Thus the function 
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[g; u, 7|[ a ,x]] by x G [a, b] is of bounded variation: 

Va[g;U,l\[a,x]\ < 00. 

18. Definition. Let / and rj be two line generalized functions on 
B(U,Av) x P(Z7) x V(U,A V ). We define a line functional with values de- 
noted by 

[lyfdv,^ 1 ®w) := [/;w\7; [t/;w,k]]| k=7 = Lfrw 1 ,*; [77; w, *]](7), 
where 7 G P(U) is a rectifiable path in U, uj^uj 1 G 25(C7, Ay) are any test 
functions. The functional J fdrj is called the non-commutative line integral 
over the Cayley-Dickson algebra A v of line generalized functions / by rj. 
Quite analogously such integral is defined for line generalized functions / 
and r] on V(U, A v ) x P{U) x V{U,A V ). 

19. Theorem. Let F and H 6e iu>o generalized functions on U , F, 5 G 
B'(U,A V ) or F, H G T>'(U,A V ), then the the non- commutative line integral 
over the Cayley-Dickson algebra A v of line generalized functions f by £ exists, 
where f is induced by F and ^ by S. 

Proof. At first it easy to mention that Definition 18 is justified by 
Definition 16 and Lemma 17.2, since the function [rj; u, k | [ 0>a; ]] is of bounded 
variation by the variable x for each rectifiable path k G P{U) and any test 
function u (see Properties 16 (Gl — G3)), while the operator / always exists 
in the class of generalized line operators, / = dg/dz, (dg(z)/dz).l = f(z) 
(see Property 16(G5)). 

Each generalized function / G B(U, A v ) can be written in the form: 

(1) [/>) = T,f^ [fj,k,^k)ij, 
where each fj^ is a real valued generalized function, fj^ G B'(U,TL), u = 
^k^k^k, G B(U, R) is a real valued test function, [fj ; k,^k) — [fj^k^k), 
[f,co) = ^2j[fj,u)ij, [fj,u)) G R for each j = 0, ...,2 V — 1 and u G B(U,A V ), 
io, ...,i2«-i is the standard base of generators of the Cayley-Dickson algebra 
A v . It is well-known that in the space B'(U, R) of generalized functions the 
space B(U, R) of test functions is everywhere dense (see [6] and §9 above). 
In view of the decomposition given by Formula (1) we get that B(U,A V ) is 
everywhere dense in B'(U,A V ). Thus sequences of test functions F l and S l 
exist converging to F and 5 correspondingly. 

Without loss of generality we can embed U into A v taking its e-enlargement 
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(open neighborhood) in case of necessity. So it is sufficient to treat the case 
of a domain U in A v . In view of the analog of the Stone- Weierstrass theorem 
(see [22j [23]) in C°(Q,A V ) for each compact canonical closed subset in A v 
the family of all super-differentiable on Q functions is dense, consequently, 
the space 14(17, Av) of all super-differentiable functions on U is everywhere 
dense in T>(U,A V ). For each rectifiable path 7 in the domain U a compact 
canonical closed domain Q exists Q C U so that 7([a, b}) C Q. Therefore, it 
is sufficient to consider test functions with compact supports in Q. Thus we 
take super-differentiable functions F n and S™. 

Let 7' be a sequence of rectifiable paths continuously differentiable, 7' G 
C x (\a, b), Av), converging to 7 in P(U) relative to the metric p. 

Then for any super-differentiable functions p and q we have 

(2) J jlP (z)dq(z) = ti(d((z)/dz).[(dq(z)/dz).d 1 l (t)}\ z=ll{t) 

= /. 6 5^ 1 (^W/^OE^ 1 («9*W/^i)dT}(*)]. 
since each super-differentiable function is Frechet differentiable, c/'yj- (t) = 

j l j'(t)dt, where (d((z)/dz).l = p(z) and for the corresponding phrases of 

them for each z G U. On the other hand, the functional 

(3) £ Et= \daz)/dz k )[T^ V = - 1 (dq k (z)/dz J )d 1 l 1 (t)] is continuous on B(U, A v ) 2 x 
P(U), i.e. for Q,p G B(U, A v ) and 7 G P(£7) as well. 

For a rectifiable path 7 in [/ it is possible to take a sequence of open e 
neighborhoods T € := \Jze-y([a,b]) B(A V , z, e), e = e(/) = 1//, where -B(v4.„, ^, e) : = 
{y : y G *4.„; |y — z| < e}. Therefore, for each function v of bounded vari- 
ation on U and each rectifiable path 7 in U a sequence of test functions 
with supports contained in T 1 ^ 1 exists such that 

]im l J u [(dC(z)/dz).e l (z)]X(dz) = I lP {z)dv{z) 
for each super-differentiable test functions]?, ( G %{U, A v ) with (d((z)/dz).l = 
p(z) on U, where A denotes the Lebesgue measure on U induced by the 
Lebesgue measure on the real shadow R 2 " of the Cayley-Dickson algebra A v , 
where 1-L(U,A V ) denotes the family of all super-differentiable functions on 
the domain U with values in the Cayley-Dickson algebra A v . 

Using the latter property and in accordance with Formulas (1 — 3) and 
16(G6, G7) we put: 

(4) [f;w,7] := limbs'; u, 7] = lim J 7 E l (y)uj(y)dy and 

44 



(5) [/V,7;H = \im[dG l /dz;co\ r ^} = lim, J 1 {(dG l (z)/dz).uj 1 (z)}du(z) 
for any 1/ G V(C/, A), where (dG l /dz).l = F l (z) on [/. 

Therefore E l converges to £ and dG l /dz converges to /, where [£; cu, *](/c|[ a ,x]) — 
[£; u), «|[ a ,x]] f° r each k G P(i7), a < x < b (see Lemma 17.2). Therefore, from 
Formulas (2 — 5) and Lemmas 17.1 and 17.2 we infer that 

(6) [J 7 /tfy, a; 1 ® u) = lim^C/ds; a; 1 , *; [S<; w , *]](V) 
= lim; jy [dG l /dz; u 1 , *; d[S'; w, *] (*)], 

where z = 7 J (t), a < t < b. 

19.1. Corollary. If F : U ^ A v is a continuous function on U and S 
zs a generalized function on U , then the non- commutative line integral over 
the Cayley-Dickson algebra A v of line generalized functions f by £ 

(i) Lf 7 M, w 1 ® w) 

exists, where f is induced by F and £ by S. 

Proof. This follows from Theorem 19 and the fact that each continuous 
function F on U gives the corresponding regular line operator valued gener- 
alized function on the space of test functions oo l in B(U,A V ) or V(U,A V ): 

[F;u\ 1 } = J J (F(z).uj 1 (z))dz. 
In this case one can take the marked function u 1 = xv, where V is a com- 
pact canonical closed sub-domain in U, since 7 ([a, b]) is compact for each 
rectifiable path 7 in U so that 7([a, &]) C V for the corresponding compact 
sub-domain V. This gives F.xv(z) — F(z) for each z £ V and F.%y(z) = 
for each z G U \ V. 

19.2. Corollary. If F e B'(U,A V ) or F G V'(U,A V ) is a general- 
ized function on U and E is a function of bounded variation on U , then the 
non- commutative line integral over the Cayley-Dickson algebra A v of line 
generalized functions f by £ 

(1) [f 7 fd£, u 1 <g> w) 
exists, where f is induced by F and £ by S. 

Proof. In this case we put 

[&w,k] := J K uj(z)dZ(z) 
for each test function oj and each rectifiable path k in U . It is sufficient to 
take marked test function oo(z) = 1 for each z G U that gives 1, *] = cG. 
Thus this corollary follows from Theorem 19. 
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19.3. Corollary. If F is a continuous function on U andS is a function 
of bounded variation on U, then the non- commutative line integral over the 
Cayley-Dickson algebra A v of line generalized functions f by £ 

(1) [J 7 fd£, u 1 ® u) 

exists, where f is induced by F and £ by 5. Moreover, this integral coincides 
with the non- commutative line integral from §5 for the unit test functions 
cu(z) = w x (z) = 1 for each z G U: 

(2) [/ 7 M,i®i) = / 7 /de- 

Proof. This follows from the combination of two preceding corollaries, 
since for a rectifiable path 7 its image in U is contained in a compact sub- 
domain V in U, i.e. 7([a, 6]) C V. 

19.4. Convolution formula for solutions of partial differential 
equations. 

Using convolutions of generalized functions a solution of the equation 

(CI) (X s +p)f = g in B(R n } Y) or in the space B'(R n , Y\ is: 

(C2) f = £rs +P *g, 
where £r a +j3 denotes a fundamental solution of the equation 

(C3) (T s + /3)£ T+/3 = 5, 
(5,0) = 0(0) (see §9). The fundamental solution of the equation 

(C4) A V = 5 with A = (X s + /?)(Tf + 
can be written as the convolution 

(C5) V =: VU = ^T«+/3 * ^x^+/3i- 

In view of Formulas 4(7 — 9) each generalized function Srs+p can also be 
found from the elliptic partial differential equation 

(C6) Sg\l/T s +/3 = ^ by the formula: 

(C7) £ r s +p = [(X s + /3)*]^t s+ /3, where 

(C8) 5^ := (X s + /3)(X s + /3)* 
(see §33 [28]). 

20. Poly-functionals. Let a fc : B(L7, A) fc -> A or a fc : £>(£/, A) fe -> A- 
be a continuous mapping satisfying the following three conditions: 
(PI) [afc, us 1 ® ... ® u k ) is R homogeneous 
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[a fc , u 1 (g...(g(u l t)(g...(guj k ) = [a fe , u) 1 <g>...<g>uj l <g)...<g)UJ k )t = [a. k t, u 1 (g...(gu k ) 
for each t G R and ^4 r additive 

[a fc , w 1 ® ... ® (a/ + «;') ® ... ® = [a fe , u 1 <8> ... ® u/ <g> ... ® + [a fc , w 1 ® 
... ® k 1 <g) ... ® 

by any A r valued test functions u l and k 1 , when other are marked, I = 1, k, 
i.e. it is k R linear and k A r additive, where [a^w 1 <g> ... <g> u k ) denotes a 
value of afc on given test A r valued functions u 1 , ...,u k ; 

(P2) [afca,^ 1 ®...® (u l f3)®...®u k ) = ([a*, co 1 ® ... <g> w' <g> ... <g> = 
[(a fc a)/3, w 1 ®...®^'®...®^) for all real- valued test functions and a, /3 G *4. r ; 

(P3) [a*,^ 1 ) <g> ... <g> = [a^w 1 <g> ... <g> for all real-valued test 

functions and each transposition a, i.e. bijective surjective mapping o : 
{1, A;} ..,*}. 

Then a^ will be called the symmetric k R linear k A r additive continu- 
ous functional, 1 < k G Z. The family of all such symmetric functionals is 
denoted by B'k, s (U, A v ) or V' ks (U, A r ) correspondingly. A functional satisfy- 
ing Conditions (PI, P2) is called a continuous /c-functional over A r and their 
family is denoted by B' k (U, Ar-) or V' k (U,A r ). When a situation is outlined 
we may omit for short " continuous" or " k R linear k A v additive" . 

The sum of two A;-functionals over the Cayley-Dickson algebra A r is pre- 
scribed by the equality: 

(PA) [a* + b fc , w 1 <g> ... <g> oo k ) = [a*, cj 1 (g) ...(g) u k ) + [b fc , u 1 (g ... (g uo k ) 
for each test functions. Using Formula (PA) each /c-functional can be written 
as 

(1) [a fc , w 1 ®...®^) = [a-kflio, w 1 <g...<gw fe ) + ... + [a fei 2'--i?2'--i, u l ®...®u k ), 
where [a fcj -, u 1 ®...®^) G R is real for all real- valued test functions u 1 , u k 
and each j; i ,...,i 2 r-i denote the standard generators of the Cayley-Dickson 
algebra A r . 

The direct product a^ <g> h p of two functionals a^ and h p for the same 
space of test functions is a k + ^functional over A r given by the following 
three conditions: 

(P5) [a* <g> bp,^ 1 <g> ... ®w fc+p ) = [a^w 1 <g> ... ® io k )[b p ,uj k+1 <g> ... <g> w fe+p ) 
for any real-valued test functions u 1 , ...,u k+p ; 

(P6) if [b p ,u) k+1 ®...®uj k+p ) G R is real for any real- valued test functions, 
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then 

[(a*iVi) <g> (b p iV 2 ), u 1 <g> ... <g> w fc+p ) = ([a* <g> b p , <g> ... <g> w fc+ P)JVi)JV 2 
for any real-valued test functions w 1 , oj k+p and Cayley-Dickson numbers 
N ± ,N 2 e A r ; 

(P7) if [a*, w 1 <g> ... <g> w fc ) G R and [b p , <g> ... <g> w fe+p ) G R are real for 
any real-valued test functions, then 

[a* <g) bp, u l <g> ... ® (w'JVi) ® <S> ^ fe+p ) = [a* <g> b p , u 1 <g> ... <g> w fc+p )iVi 
for any real- valued test functions w 1 , ...,u k+p and each Cayley-Dickson num- 
ber Ni G v4 r for each I — 1, k + p. 

Therefore, we can now consider a partial differential operator of order 
u acting on a generalized function / G B'(U,A r ) or / G V'(U,A r ) and 
with generalized coefficients either a a G #| Q |(£/, „4 r ) or all a a G P'| a |(C/, Ar) 
correspondingly: 

(1) Af(x) = £ (d a f(x)) ® KM*)) ® 1 0(uHq|) ], 

|a|<u 

where <9 a / = <9' Q '/(a;)/(9a;o ...(9x" n , a; = xo^o + ■■■x n i n , Xj G R for each 
j, 1 < n = 2 r — 1, ct = (cto, ct n ), | or | = «o + ••• + < a-,- G Z, 
[1, w) := u>(y)X(dy), A denotes the Lebesgue measure on U, for convenience 
1®° means the multiplication on the unit 1 G R. The partial differential 
equation 

(2) Af = g in terms of generalized functions has a solution / means by 
the definition that 

(3) [Af,io^ u+1 ^) = [g,uj® (u +V) 

for each real- valued test function u on U, where u® k = u <g> ... <g> u denotes 
the k times direct product of a test functions u. 

21. Theorem. Let A = A u be a partial differential operator with general- 
ized over the Cayley-Dickson algebra A r coefficients of an even order u = 2s 
on U such that each a a is symmetric for \a\ — u and A has the form 

(4) Af = (B U:1 f)c Ujl + ... + (B U:k f)c Ujk , where each 
B u ,p B U p Q -\- Qu—i,p 

is a partial differential operator by variables x muA+ ... +rn ^ p _ 1+ i,...,x mul+ ... +mup 
and of the order u, m Uj0 = 0, c Ujk {x) G A r for each k, its principal part 
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(6) B u , pfi f = Ei a \=s(d 2a f)®* P ,2a(x) 

is elliptic, i.e. EH= s 2/ 2a [ap, 2a ,w® 2s ) > for all y k(1) ,...,y k(muj>) in R with 
fc(l) = m tt) i+...+m U) p_i+l,...,A;(m U) p) = m U)1 +...+m U)P , y 13 := yilly--y k [^) 
and [a Pj2Q ,, w® 2s ) > /or eac/i rea/ £est function u, either < r < 3 and f is 
with values in A r , or r > 4 and f is real-valued on real-valued test functions. 
Then three partial differential operators X s and T\ and Q of orders s and p 
with p < u — 1 with generalized on U coefficients with values in A v exist such 
that 

(7) [Af,u®( u+ V) = [T S (T S J) + Qf,uj^ u+1 ^ for each real-valued test func- 
tion u on U . 

Proof. If a 2s is a symmetric functional and [c S) us® s ) = [a 2s , a;® 25 ) 1 / 2 for 
each real-valued test function u, then by Formulas 20(P1, P2) this functional 
c s has an extension up to a continuous s-functional over the Cayley-Dickson 
algebra A r - This is sufficient for Formula (7), where only real- valued test 
functions u> are taken. 

Consider a continuous p-functional c p over A v , p G N. Supply the domain 
U with the metric induced from the corresponding Euclidean space or the 
Cayley-Dickson algebra in which U is embedded. It is possible to take a 
sequence of non-negative test functions ioo on U with a support supp(iu) 
contained in the ball B(U, z,l/l) with center z and radius l/l and iu positive 
on some open neighborhood of a marked point z in U so that f v iuu(z)X(dz) = 
1 for each I G N. If the p-functional c p is regular and realized by a continuous 
A v valued function g on U p , then liim[c p , cu® p ) = g(z, z). Thus the partial 
differential equation 20(2) for regular functionals and their derivatives implies 
the classical partial differential equation 2(1). 

Therefore, the statement of this theorem follows from Theorem 10, and 
§§14, 15 and 20, since the spaces of test functions are dense in the spaces of 
generalized functions (see §19). 

22. Corollary. If Af = E J , k (d 2 f(z)/dz k dz J )m J , k (z)+E 3 (df(z)/dz ] )^ 
bj(z) <S> 1 + f(z) <S> r)(z) <8> 1 is a second order partial differential operator 
with generalized coefficients in B'(U,A r ) or V'{U,A r ), where each a,j tk is 
symmetric, f and A r are as in §20, then three partial differential operators 
T + (3, Ti + j3i and Q of the first order with generalized coefficients with 
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values in A v for suitable v > r of the same class exist such that 

(1) [Af,u m ) = [(T + f3)({T 1 + f3 1 )f + Qf),u m ) for each real-valued test 
function u on U. 

Proof. This follows from Theorem 21 and Corollary 12 and §§2 and 8. 

23. Anti-derivatives of first order partial differential operators 
with generalized coefficients. 

Theorem. Let T be a first order partial differential operator given by the 
formula 

(1) T/ = T,]=o(df/dzj) <S> [iji/>j(z)] or 

(2) rf = E] =0 (df/dz 3 )®<p*(z), 

where supp(ipj(z)) = U or supp((f>j(z)) = U for each j respectively, f and 
ipj(z) or if)j(z) are A v -valued generalized functions inB'{U,A r ) orV(U,A r ) 
on the domainU satisfying Conditions 1(D1,D2), alg^{[(f)j, uu), [(f) k ,uj), [(pi,oj)} 
is alternative for all < j, k, I < 2 V — 1 and a/gR,{[0O) w )j •••■> [02 v -i, w )} C A v 
for each real-valued test function uj on U . Then its anti-derivative opera- 
tor Xt exists such that TX^f = f for each continuous generalized function 
f : U — >■ A v and it has an expression through line integrals of generalized 
functions. 

Proof. When an operator with generalized coefficients is given by For- 
mula (1), we shall take unknown generalized functions fj(z) G A v as solutions 
of the system of partial differential equations by real variables z^. 

(3) [(dvj(z)/dzj) <g>^(z),w® 2 ) = [l,u m ) for all 1 < j < n; 

(4) [ij, k {z) ® (du 1 (z)/dz k ),u^) = tyj(z) ® (dis k (z)/d Zj ),u® 2 ) for all 1 < 
j < k < n and and real-valued test functions u on U . 

If the operator is given by Formula (2) we consider the system of partial 
differential equations: 

(5) [{{dg{z)ldz)\dvj{z)ldz k \)®^ = 
for all < j < k < n; 

(6) dvj(z)/dzj = 4>j{z) for all j — 0, n\ 

(7) [{{dg{z)/dz)4j{z))®(t)*{z),uj® 2 ) = [/(z)®l,w® 2 ) for each j = 0,...,n 
and every real- valued test function uj. 

Certainly the system of differential equations given by Formulas (3, 4) or 
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(5 — 7) have solutions in the spaces of test functions B(U,A r ) or V(U,A r ), 
when all functions ipj or <pj are in the same space respectively. Applying 
§§4 or 5 we find generalized functions Uj resolving these system of partial 
differential equations correspondingly, when all functions ipj or <fij are gen- 
eralized functions, since the spaces of test functions are dense in the spaces 
of generalized functions (see §19). Substituting line integrals J^q(y)duj(y) 
from §§4 and 5 on line integrals [f q{y)duj{y) 1 uj 1 <8> u) from §19 one gets the 
statement of this theorem, since test functions u 1 and u in the line integrals 
of generalized functions can also be taken real- valued and the real field is the 
center of the Cayley-Dickson algebra A v . Therefore, we infer that 
(8) d[L„| f(y)dvj(y),uj®uj)/dz k = [f(z).[du j (z)/dz k ],u®uj) 

I I Aft jlz J 

for each real-valued test function u and z G U ', where j a (t z ) = z, t z G< 
a a ,b a >, a G A. Equality (8), Theorem 19 and Corollaries 19.1-19.3 and 
Conditions 20(P1 — PI) give the formula for an anti-derivative operator: 

(9) [l T f,u®u>) = [ r I f{z)dz,u®u) = (n+iy^ilf q{y)dv 3 {y),u®u) 

J j=0 J -1 a \[a a ,t\ 

for each real- valued test-function co, where a G A, a a < t < b a , t = t z , 
z = j(t), consequently, 

(10) [T T / f{y)dy, uj^) = [f ® 1 ® 1, ^ 3 ). 

23.1. Note. Certainly, the case of the partial differential operator 

(1) Tf = T$= Q (df/dz h(J) )®<l>i (j) {z), 
where < k(0) < k(l) < ... < k(n) < 2 V — 1 reduces to the considered in §23 
case by a suitable change of variables z !->■ y so that z k m = yy 

24. Example. We consider a consequence of Formulas 15(2 — 6). If q(t) 
is a differentiable function on the real field R having simple zeros q(tj) = 
(i.e. zeros of the first order), then 

(1) S(q(t)) = E J ^ ) S(t-t J ), 

where the sum is accomplished by all zeros tj of the equation q(t) =0 (see 
Formula II.2.6(J^) [6]). Therefore, if j(t) is a C 1 path in U intersecting 
the surface dU s D dU p at the marked point x of index I = l SjP (x), 7(t ) = x, 
< r < 1, such that j(t) G U s for each r < r and j(t) G U p for each 
r > t then 

(2) dg^/dr = dg'i^/dr + h£* =1 S{P j )dP j { 1 {r)) / dr, 
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where 0(t) = for t < and 6(t) = 1 for t > 0, ^(r)) = 3(7(7-)) - h9(r - 
To), 

(2.1) /i = lim r ; ro g{i{r)) - lim TtT0 0(7(7-)), 

(2.2) dPjijir^/dr = TJU(dP 3 (z)/dz k )(d lk (r)/dr)\ z=l[T) (see Example 
1.2.2.2 [6]). Particularly, if a point a; is of index 1, then Formula (2) simplifies: 

(3) dg{^{r))/dr = ^ 1 ( 7 (r))/rfr + /w5(P)[dP(7(r))/dr]. Particularly, 
these formulas can be applied to dvj. 

Let a partial differential differential operator Q be given by Formula 14(7) 
and functions v k are found (see Theorems 5 and 21 above). We put in 
accordance with Formula 15(6) 

( 4 ) v\(z) = u k( z ) - T,s, p hk;s >P {z)9{P SjP .i{z),...,P s ^ m .{z)), 
where 

(5) h k . S)P (z(x)) = h k (x)\ xe rs,p := 

'in'v::- ,.!' «....,.o v k (z(x + y)) - liin, /:; _ , .o u k (z(x - y)), 

where the sum is by s and p with dU s D dU p 7^ 0. 

Let a domain be a canonical closed compact set in the Euclidean space 
R n+1 embedded into A v and contained in a canonical closed compact domain 
U so that W = {z G U : z,- = \/n < j < T - 1}. Thus T from test and 
generalized functions on W is extended on test and generalized functions on 
U. We can put Uj = for n < j < 2 V — 1, when n < 2 V — 1. Then for the 
rectifiable path 7 (see above) we get 

(6) (n + ir 1 ^/ q(y)du k (y),u 1 xu) 

k=0 J Tl[a,t] 

= (» + l)- l £{[/ q(y)dvl(yW x u) 

l n 

+ [<?;^ 1 ,*; EW^E^-Pj) H(^W/^m)(57m(7-)/rfr)| z=7(r) ;u;,*]](7)}, 

s,p j=l m=l 

where 7 G {7° : a G A} is taken from the foliation C 1 family of paths 
(see §6.1 above and also Theorem 2.13 [27]), z = 7(^2), t z G [a a ,b a ], [a,t] = 
[a a , t z ], I — l s ,p(z) denotes an index of a point z in the intersection of bound- 
aries T s ' p := dU s fl dU p 7^ 0, and u are real- valued test functions. Since 
uj 1 is real- valued, we get /(z).^ 1 = f(z)cu 1 (z) and 

(7) [q; u>\ *; [Es, P hk;s, P (z) Ej=i WO Em=i ( 9P j(z)/ 9z m ) (dj m (r) / dr ) | *= 7 ( T ) ; w, *]]( 7 ) 
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= E^fe^u; 1 ^); [0(t - r SjP )h k . BtP (z),uj))\ z=y{T) , 
where t SjP corresponds to the intersection point 7(t SiP ) of 7 with r s,p 7^ 0. 
Here the expression [q, u)\ z=1 ( T ) := lining o k^oj o «;■?) denotes the restriction 
of the generalized function from U onto 7([a, 6]), k;- 7 G X>(i7, Ay) is a sequence 
of test functions and K?((f>([a,b])) C U for each j e N, G X>([a, &], A„), 
Djli supp(K j ) = 0([a, 6]), linx, k- 7 o = 7 in ~P(U). Therefore, the derivative 
of the operator [(n + l) -1 Efc=o[/ 7 | [ot] Q^dvkiy),^ 1 x 00) by the parameter 
r G [a, b] for the real test functions u 1 and w is the following: 

(8) d(n+ l^ELoUyi^ q(y)^ k (vW x w)/0r = 

[(n+l)- 1 £jU{l®?H*)-(*'J^^ 

W(7(r))/rfr))| z=7(r) ) <g> 5(r - r s , p ), <g> w), 

where dg(z)/dz = q(z) on [/ in the class of generalized operator L ? (A„, A„) 
valued functions, (dg(z)/dz).l = q{z) on [/, h g s \ p (z) = h(z) is given by For- 
mula (2.1) for the derivative operator dg(z)/dz = g' instead of g on each 
T s ' p 7^ 0, q 1 is given by Formula 15(6) for the function q(z) with values in 
L q (A v: A v ) instead of g(z). The terms like q 1 (z).(dul('y(T))/dT) correspond 
to the action of the operator valued generalized function ^(z) on the gener- 
alized function (dv>l(j(r)) / 'dr) which gives a generalized function. 

Using Formulas (6 — 8) for n constant on U and ipj(z) or (f)j(z) respectively 
non-zero for each z G U and all j — 0, n we infer that for a continuous or 
generalized function / 

(9) Tlrf(z) = f(z), where 

(10) T / f(z)dz := {(n+1)- 1 J2 I q(z)duj(z), a G A, a a < t < b a }, 

J j=0 J ~f a \[a a ,t] 

where q = (dg/dz).l and g is given by the Equation 5(9), since f 1 = f and 
h 9 s ' p = in the class of generalized functions / and in the class of continuous 
functions /, also hk- s ,p = for = v\ on U. 

Formulas (9, 10) show what sort of boundary conditions is sufficient to 
specify a unique solution for a given domain U with sub-domains U s . If U is 
C 1 diffeomorphic to the half-space H p := {z G A v : ZoPo + ... + z 2 > 
0}, where p = p i + ... + p 2 "-iZ2"-i is a marked Cayley-Dickson number, 
Po, ...,p 2 v_i G R, and sub-domains U s are not prescribed, then it is sufficient 
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to give the boundary condition F\qu = G when a solution is in the class 
of continuous or generalized functions with the corresponding / and ipj or 
<f>j. Indeed, if the functions v k along 7" are defined up to constants /ik, 
the differentials are the same d[yk + IJ>k){ z )\z=-y a (T) — dvk\ z =^ a (r) i n the anti- 
derivative operator, when rf/Xfc| 2=7 c«( T ) = for each a G A and r G [a a , b a ]. 

The operator Xt may be applied also piecewise on each U s . If a solution 
F is locally continuous on U and continuous on each sub-domain U s , then 
boundary conditions F\g Ua = G s for all s = l,...,m may be necessary to 
specify a solution F. Without boundary conditions the anti-derivative oper- 
ator applied to / gives the general solution lyf of the differential equation 
TF = f. 

If each u k is continuously differentiable, which is possible, when each 
function ipk or (pk is continuous, and / is continuous on U, then a solution 
F — {t J 7 «| [o t] f(y)dy '■ a E A, a a < t < b a } is continuously differentiable 
by each z k , z e U, z = 7 a (t 2 ). 

One can also mention that the sequence m uj(zo : z n ) = (27rm) _ (™ +1 ^ 2 exp{ — (z^+ 
... + z^)/(2m)} converges to the delta-function on R n+1 embedded into A v , 
m e N. Then the sequence m 6(z ) = Jl° m (27rm)~ 1/2 exp{-t 2 / (2m)}dt con- 
verges to the function, while each function m 6(zo) is analytic by Zq, since 
m G N, exp{— t 2 /(2m)} is the analytic function with the infinite radius of 
convergence of its power series, while lim m ^ +00 /J™(27rm) -1 / 2 exp{— t 2 / (2m)}dt = 
0. Then each m u(zo, z n ) and l\j=o m@( z j) can be written in the ^-representation 
over A v as the analytic function with the help of Formulas 1(1—3), where n < 
2 V —1, zj G Rfor each j, z = z i +...+z 2 ^-ii2^i- Thus limm^oo m u(zo, z n ) = 
5(z , ...,z n ) and lim^^ n" =0 mO(zj) = U]=o°( z j)i d9(z )/dz = 5(z ). 

25. Boundary conditions. 

If U is a domain as in §15, we put B(dU, Y) = {f\ 9u : / G B(U, Y)} and 
T>(dU, Y) — {f\ou '■ f £ T^(U, Y)} when a boundary dU is non-void so that 
the topologically adjoint linear over R spaces, left and right A r modules, of 
generalized functions are B'{dU,Y) and V'(dU,Y). 

Let us consider a generalized function / G B'(dU, Y) or V'(dU, Y) and 
a test function h G B(dU, Y) or V(dU, Y) respectively. One can take g G 
B(U, Y) or V{U, Y) and a sequence q m G B(C/, F) or V(U, Y) with supports 
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non intersecting with the boundary supp(q m ) fl dU = such that (g — q m ) 
tends to zero in B(V, Y) or V(V, Y) for each compact subset V in the interior 
Int(U), when m tends to the infinity, while lim m (<? — q m ) = h in B(dU,Y) 
or T>(dU, Y) respectively Here as usually the interior Int(U) is taken in the 
corresponding topological space R ra or A r . Each generalized function is a 
limit of test functions, consequently, a generalized function £ G B'(U,Y) or 
£>'([/, y) exists so that 

(51) Km m [Z,(g-q m )) = [f,h). 
Vise versa if £ G y) or T>'(U, Y) is a generalized function on U, then For- 
mula (Bl) defines a generalized function / G B'(dU,Y) or V'(dU,Y), which 
we consider as the restriction of £ on B(dU, Y) or V(dU, Y) correspondingly. 
In view of the definition of convergence of test and generalized functions For- 
mula (Bl) defines the unique restriction / for the given generalized function 

A subsequent use of decomposition of operators into compositions of first 
order partial differential operators and their anti-derivation operators permits 
to integrate partial differential equations with locally continuous or general- 
ized coefficients. 

The results and definitions of previous sections show that for the differ- 
ential equation 

(1) Af = g, 

where a partial differential operator is written in accordance with Formulas 
10(1,2). When dU is a C 1 -manifold without corner points of index greater 
than one, the following boundary conditions may be used: 

(2) f(t)\au = hit'), (dMf(t)/dsf...ds%)\ du = f {q) (t>) for \q\ < a - 1, 
where s = (s u s n ) G R n , (q) = (q u q n ), \q\ = q 1 + ... + q n , < q k G Z 
for each k, t G dU is denoted by t' , /o, f( q ) are given functions. Generally 
these conditions may be excessive, so one uses some of them or their linear 
combinations (see (4) below). Frequently, the boundary conditions 

(3) f(t)\ du = f (f), (d l f(t)/dv l )\ 9u = hit') for 1 <l< a - 1 are also 
used, where v denotes a real variable along a unit external normal to the 
boundary dU at a point t' G dll . Using partial differentiation in local 
coordinates on dU and (3) one can calculate in principle all other boundary 
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conditions in (2) almost everywhere on dU . 

It is possible to describe as an example a particular class of domains and 
boundary conditions. Suppose that a domain U\ and its boundary dU\ satisfy 
Conditions (Dl,i — vii) and g 1 = gxui is a regular or generalized function 
on R n with its support in Ui. Then any function g on R n gives the function 
9Xu 2 ='■ 92 on R™, where U 2 = R n \ U\. Take now new domain U satisfying 
Conditions (Dl,i — vii) and (D2 — D5): 

(D2) UDUr and dU C dUi, 

(D3) if a straight line £ containing a point W\ (see lb(vi)) intersects dU\ 
at two points y\ and 1/2, then only one point either y\ or y 2 belongs to dU, 
where w\ G U±, U — w\ and Ui — w\ are convex; if £ intersects dU\ only at 
one point, then it intersects dU at the same point; 

(DA) any straight line £ through the point w\ either does not intersect 
dU or intersects the boundary dU only at one point. 

Take now g with supp(g) C U, then supp(gxui) C C/i. Therefore, any 
problem (1) on C/i can be considered as the restriction of the problem (1) 
defined on U, satisfying (Dl — DA,i — vii). Any solution / of (1) on U with 
the boundary conditions on dU gives the solution as the restriction f\ Vl on 
U\ with the boundary conditions on dU . 

Henceforward, we suppose that the domain U satisfies Conditions (Dl, DA, i— 
vii), which are rather mild and natural. 

Thus the sufficient boundary conditions are: 

(4) (^i^/ctw))/^^...^^)!^^ = ^c)^)) 

for \(3\ = \q\, where m = \h(lj)\, \j\ < a, \(lj)\ > 1, ^ 0, q k = for 
hjk = 0, m k + q k + h k = j k , h k = sign(l k j k ), < q k < j k - 1 for k > n - re; 
4 > /3,(i)(t^) are known functions on dU^, G dUyy In the half-space t n > 
only the partial derivatives by t n 

(5) d^f(t)/dt\ tn=0 

are necessary for j3 — \q\ < a and q as above. 

Depending on coefficients of the operator A and the domain U some 
boundary conditions may be dropped, when the corresponding terms vanish. 

Conditions in (5) are given on disjoint for different (/) sub-manifolds dU^ 
in dU and partial derivatives are along orthogonal to them coordinates in R n , 
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so they are correctly posed. 

We recall, that a characteristic surface of a partial differential operator 
given by Formula 10(1) is a surface defined as a zero of C u differentiable 
function (f)(xi, x n ) = in the Euclidean space or in the Cayley-Dickson 
algebra such that at each point x of it the condition is satisfied 

(cs) T.\ 3 \=u^Mx)){d<p/dx 1 y\..(d(p/dx n y- = o 

and at least one of the partial derivatives {dip/ dxu) 7^ is non-zero. Generally 
a domain U is worthwhile to choose with its interior Int(U) non-intersecting 
with a characteristic surface 0(xi, ■ ■■,x n ) = (see also [3"2"| 157]). 

26. Solutions of second order partial differential equations with 
the help of the line integration over the Cayley-Dickson algebras. 

Mention that the operator (T+/3)(z , z n ) acting on a function depend- 
ing on variables zq, ...,z n only can be written as 

(1) *(*b, z n+1 )(f(z)z n+1 )\ Zn+1=1 = T(z , z n )f(z) + f(z)P(z) 
= [Y,%l{d{f{z)z n ^)/dz j )^ j {z)\\ Zn+1=1 , 

where (fi* n+1 (z) = f3(z), each function <f>j(z) and f(z) may depend on z , z n 
only, omitting for short the direct product (g) in the case of generalized coeffi- 
cients in Formula (1) and henceforth. The operator \l/(zo; z n+i){f{z)z n +i) 
may be reduced to the form satisfying conditions of Theorems 5 or 23 using a 
suitable change of variables. This procedure gives an anti-derivative operator 

(1.1) X T+/ 3 =Z*L +1 =i such that 

(1.2) (T + p)Z r+ pf = f 

for a continuous function or a generalized function /. Therefore, we shall 
consider operators of the form T and their compositions and sums. 

We take the partial differential equation with piecewise continuous or 
generalized coefficients 

(2) A = TTi/(z) + T 2 f(z) = g, 
where 

(3) T k f{z) = \&aPf{*)ld*iW*Y\ 

for k = 1,2 or without this index k, <p k j{z) = ijipj{z) (see §§4, 5 and 23). 
For solving it we write the system: 

(4) T 1 f = g 1 ,Tg l = g-T 2 f. 

In accordance with Equation (2) we have: 
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(5) Tg 1 + T 2 (Ti 1 g 1 )=g, 

where the inverse operator Tf 1 is the ant i- derivation operator X Tl described 
above in Theorems 4, 5 and 23. If T 2 ^ we suppose that either (i) gi or 
[gi,uj) is real- valued or the algebra alga.{gi(z), 4>}(z), <frl(z)} or aZ# R {[g>i, u), [0}, w), w)} 
for all j, = 0, ...,n is associative for each z G U in the case of functions or 
for every real- valued test function a; in the case of generalized functions cor- 
respondingly. Calculating the expression 

(6) Tg 1 + T 2 (l ri g 1 ) = (T + ^)g 1 = g 

we get a new operator (T + /3 3 ) = ^ in accordance with Formulas 4(5,6), 
5(6 — 8) and 23(3 — 7) omitting (g) and ®u and [*, u® 3 ) for short in the case 
of generalized coefficients, where 

when a solution {vj(z) : j} is chosen real, i.e. each function Vj is real- valued 
or a real- valued generalized function on real- valued test functions (see System 
4(10,11) and §23). 

Generally without supposition {%) we deduce that 

(6.1) Tg hl + T 2 (X Tl5M ) = (T + = itf 

for each / = 0, ...,n, where f3f = «*[E",jfc=o(*j(^j /^fc))(0l)1> 

(6-2) JZo!? = 9> 

(6.3) £p=o0i,J*J = #i, 

each <7i,j(z) or [<7i,j,u;) is real-valued for each z G £/ or every real-valued test 
function oj respectively. Solving the system (6.1 — 3) with the help of known 
anti-derivative operators one finds gi, where VPj = T + j3f. Thus an 
anti-derivative operator J : = 1 7t+t 2 x Ti exists so that 

(6.4) (T + T 2 X Tl ) Jr+rn ri 9 = 9- 

In the particular case (i) the equality J7t+t 2 x Ti = ^* is satisfied. 

Therefore, in the case of either continuous coefficients of operators and g 
or generalized coefficients and g the general solution is: 

(7) / = X Tl gi = 1 Tl Jg, where 

(8) gi = Jg- 

If T 2 = the formula simplifies to 

(9) / = l Tl l r g. 

Examples of boundary conditions and domains permitting to specify a 
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unique solution are given in §25 above. 

27. Example. Let us consider a function and its phrase satisfying 
Condition 7(P3). Therefore, we get 

(1) T^( 9l )(z a i(x)) = -<j> f ,(Im(z a i(x))) + f: gi (z)dz + <f> gi (x>) 
in accordance with Formula 7(6), where x' and Im(z a i(x)) are written in the 
x and ^-representations respectively using Identities 1(1 — 3). In particular, 
for T 1 = <j\ we have the coefficients if>l(z) = ii(sL k (z)wjilr k ) for each / = 2 r k 
with k G {mi + ••• +m,j-i + 1, m 1 + ... +m,j}, while ipf(z) = for all others 
I for each z. A function a 1 is given by Formula 7(3) for if) 1 instead of if). Let 
the first order operator Q be written in its standard form: 

m m 1 +...+mj l 

( 2 ) Qf = 12 12 s k (z)(df/dz 2 r k ){u j i 2 r k ), 

j=l k=mi+...+mj l -i+l 

since iykWj = w*%i-r k for each Wj G A r and k > 1, when r > 1, where 
iUj G »4 r and -u^ = tij(-z) G A„ for each j, s k (z) is the real- valued (super- 
)differentiable function for each k. If £ = 2 (y), then (d£ / dy) .[(dy / d£) .h] = h 
for each Cayley-Dickson number ft, G Au- This implies that these two R-linear 
A„-additive operators are related by the equality (dy/d£).h = 'dy)^ 1 .h. 
On the other hand, (dz a (y) / dy) .1 = a(y) G A^ and y E U C A v in the 
considered case. The function z a (y) is defined up to the addendum z a (a), 
where a G H ao fl £/, if ao := {z G Ay : ite(z) = Re(a )}. We can choose 
a (y') s° that (dx/dy).(l/a) = 1 for each y for which a = a(y) ^ and 
inevitably we get (dy/dx).l = l/a(y). 

In the particular case of a, o\ and Q accomplishing the differentiation 
with the help of the latter identities we infer that: 

m m 1 +...+m jl 

(3) Qfc 1 gi)(x) = -12 12 [{d<t> r {x')/dx).ir k 

j=l k=mi+...+m,j 1 -i+l 

+{9i(y) + (y')/dy)}.[(dz a i (y) / dy)' 1 .i 2 r k \ (s k (y)uj(y)i 2 r k ), 

where g\ = d(i(y)/dy for a (super-) differentiable function d such that (d(±(y) / dy) .1 = 
gi(y'), i>l(z) = ii(sL k (z)wji2r k ) for each / = 2 r A; with G {mi + ... + m.,_i + 
1, mi + ... + rrij}, while if>j(z) = for all others / for each z. Also ipi(z) = 
i[(eL k (z)w*i2r k ) for each / = Tk with k G {TOi + ...+TOj_i + l, ...,m 1 + ...+mj}, 
while if>i(z) = for all others / and for each z. We introduce the notation: 
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(4) 6i(y) = i l [(dz a (y)/dy) 1 .i 2 r k }{s k {y)u j (y)i 2 r k ) for / = 2 r k and k G 
{mi + ••• + m ji-i + 1, + ••• + m^}, ^(?/) = f° r a ^ others I and for 

each y; 

m rm+.-.+mj 

(5) a(x) = —y~]w.j afc(x) = cr(x) and 

jr'=l fe=mi+...+m J -_i+l 

(6) «j(x) = ii[sk(x)uj(x)i2rk] for / = Tk with fc e {mi + ... + m 3x -\ + 
1, mi + ... + mjj}, Ki(x) = for every other I and for each x. 
In the general case 

(?) r 2 (T^ gi )(x) = -i;m f ,v)/dx).i j 

+{9i(y) + (d^iy^/dy^.lidZa^/dyyKij^jix)), 

where g\ = d(i(y)/dy for a (super-)differentiable function d such that (d(i(y) / dy) .1 = 
9i(y')- We shall use the notation: 

(8) 9j(y) = i j [(dz a (y)/dy)- 1 .i j ](i*ip j (y)) and for each y and each j; 

(9) a(x) = T(x) and 

(10) «j(x) = ij[ijfpj(x)] for each x and j. Substituting (7) into 26(5) we 
deduce that: 

(11) gi(z a (x)) +g 1 (z e (x)) = ((x), where 

rx 

(12) ((x) = -(f)g 1 (Irn(z a (x)))-(j)g 1 (Im(z0(x)))+(j)ff(Im(z K (x)))+ g(z)dz. 

J a 

For x = a we certainly have z a (a) = ze(a). Suppose that a(x) ^ 6{x) 
identically. The dimension of the Cayley-Dickson algebra A v over the real 
field is not less, than four. Therefore, we can choose a path 7 so that 7 
is orthogonal to 9 and k at each point on 7, that is "/(t) _L 6'{j{i)) and 
■y'(t) _L /t'(7(t)) relative to the real- valued scalar product (RS) for each t e 
[0, 1], where ^'(t) := d^(t)/dt. Then <?i(z0(x)) = <?i(a) and (j)f>(Im(z K (x))) = 
4>f(a') for each x = ^(t). Therefore, along such path 7 one has 

(13) g l (z a (x))+g l (a) = C(x) = -(f) gi (Im(z a (x))-(f)g 1 (a')+<f) f (a')+ f g(z)dz 

Ja 

for each x = "f(t). Expressing g±(z) from Equation (11), substituting into 
26(4) and integrating one gets: 

(14) f(z a i(x)) = -<t) f (Im(z a i(x))) + / gi (z)dz. 

Ja 
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Particularly, if the operator A is with constant coefficients, then s^(x) = 
identically for each k, consequently, 8 = and k = identically and 
— 9i{ a ) — 4>gi{I m {z a {x))) for each x, when / has a right linear 
derivative by z. Arbitrary integration terms in (11, 14) can be specified from 
the boundary conditions. 

Finally, the restriction from the domain in A v onto the initial domain of 
real variables in the real shadow and the extraction of K v r o / g A r with the 
help of Formulas 1(1 — 3) gives the reduction of a solution from A v to A r , 
where tt£ : Av — > A r is the R-linear projection operator defined as the sum 
of projection operators ttq + ... + ^2 r -i given by Formulas 3(P1, P2) on Ri,- 
for j = 0,...,2 r - 1. 

28. Laplace's operator. When 

(1) A = A n = E] =1 d 2 /dzl 
is Laplace's operator, then 

(2) Tf(z) = ZU(df(z)/dzj)i*, so that 

(3) A n = TT* = -TT, T 1 = -T, 

where 2 < n < 2 r — 1, zi, ...,z n G R, in accordance with §2. Consider the 
fundamental solution \l/ of the following equation 

(4) E^(zi, z n ) = 5(zi, z n ) with 5 = A n satisfies the identity: 

(5) \1/ = — (T\&) * (T\I>) (see the convolution of generalized functions and 
this formula in §§9, 19.4). 

We seek the real fundamental solution \1/ = ty n , since the Laplace operator 
is real. The Fourier transform with the generator i (see §33 [28]) by real 
variables Z\, z n gives 

(6) F(* n )(x) = -[F(T^ n )(x)] 2 = -l-i(E]=ixrt)F(V n )(x)] 2 , since 
F(T* n )(x) = TJ; =l F{d^ n /dz 3 )i* = -i(ZU x ^) F ^n)(x), 

where x = (x±, ...,x n ), x±, ...,x n G R (see also §2). Thus we get the identity 

(7) F(* n )(x) = -{Y.U x W^n){x)f or 

(8) F(\E' n )(x) = —(1/(^=1^)) for n > 3 is the regular generalized 
function (functional), while 

(9) F(* 2 )(x) = -P(l/(E?=i^ 2 )) for n = 2. 

We recall that the generalized function V(l/(J2 2 =i x2 )) on G X?(R 2 ,R) is 
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defined as the regularization: 

2 

(V(l/(£^)),<f>) = I [<P{x)-m\\A- 2 dx+ I <P{x)\x\~ 2 dx, 
~i A*\<i M>i 

where x = (xi,x 2 ), \x\ 2 = x\ + x 2 , X\,x 2 G R. 

The inverse Fourier transform (F~ 1 g)(x) = (27i)~ n (Fg)(—x) of the func- 
tions l/(E"=i z]) for n > 3 and P(l/(£j=i 2?)) for n = 2 in the class of the 
generalized functions is known (see [6] and §§9.7 and 11.8 [37]) and gives 

(10) * n ( Zl ,...,z n ) = Cnfe^z]) 1 -^ 2 for 3 < n, where C n = -l/[(n - 
2)a n ], a n = 4n n / 2 /T((n/2) — 1) denotes the surface of the unit sphere in R n , 
T(x) denotes Euler's gamma-function, while 

(11) * 2 ( Zl , z 2 ) = C 2 ln(E?=i zf) for n = 2, where C 2 = l/(4vr). 

Thus the technique of convolutions over the Cayley-Dickson algebra has per- 
mitted to get the solution of the Laplace operator. 

Another method is with the line integration over the Cayley-Dickson al- 
gebras. In accordance with Formula 26(9) we get 

ty n (z 1 ,...,z n ) = -l r l r 5. 

Laplace's operator and the delta-function are invariant under any orthog- 
onal transform T G 0„(R) of R n . Therefore, a fundamental solution also 
is invariant relative to the orthogonal group 0„(R). That is \l/ n depends 
on \z\ and is independent of spherical angles in spherical system of coor- 
dinates. Thus we choose the corresponding branches of the anti-derivative 
ZyZy<5. The volume element in the Euclidean space R n can be written as 
X(dz) = x n ~ l dxds ) where x = \z\ and ds is the surface element (measure) on 
the unit sphere S*™ -1 . For each orthogonal transform its Jacobian is unit. 

One can take the family of test functions r/ e = ( 27r )i/2 e n ex P{~( z i + ••• + 
z 2 )/(2e 2 )} tending to the delta-function, when e > tends to zero. These 
functions can be written in the z- representation due to Formulas 1(1 — 3). On 
the other hand, for each z-analytic function rj with real expansion coefficients 
into a power series each line integral over the Cayley-Dickson algebra A v 
restricted on any complex plane R © MR coincides with the usual complex 
integral, where M is a purely imaginary Cayley-Dickson number. Therefore, 

JlU!v(z)dz]y k dy = ^ J^ 1 - z k ^) V (z)dz for k ? -1 and 

Sl\J2v(*)dz]$dy = tiiHt) - Hz)) V (z)dz. 
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For the characteristic function XB(R n ,o,x) of the ball -B(R n ,0,:r) of radius 
x > with the center at zero in the Euclidean space R n embedded into 
the Cayley-Dickson algebra A v one can take a sequence of test functions 
iu 1 converging to the regular generalized function Xb(R",o,i)) when / tends 
to the infinity, consequently, linn J R „ iUJ 1 (z)X(dz) = a n x n . Integrating twice 
with the anti-derivative operator these test functions i] e in accordance with 
Example 4.1 and Theorems 19 and 23 and taking the limit with e tending to 
zero from the right one gets Formulas (10, 11). 

This can also be deduced with the help of the Fourier transform with the 
generator i: 

(12) F{* n )(x) = F(-l r l r 6) = {YT j =ixf)- 1 F(5) = (E"=i^) -1 - 
Applying the inverse Fourier transform to both sides of Equation (12) we 
again get Formulas (10, 11). 

29. The hyperbolic operators with constant coefficients. 

Consider now the hyperbolic operator 

(1) A = L M = EU - £i= P+ i &/9*h 

where p + q = n, 1 < p and 1 < q, (p, q) is the signature of this operator, 
Z\,...,z n G R. Take two operators T and T 1 with constant A v coefficients so 
that 

(2) Tf(z) = EUidfW/^h + ^]= P+ i(9f(z)/dz j )[iti* 2j ] and 

(3) T!/(z) = n = i(9/(«)/^)<i- + E^.i(9/(«)/^)[<i<y, so that 

(4) L M = TTi, 

where 2 < n < 2 v ~ r — 1, r — 1 < v, in accordance with Formulas 2(7 — 9). 
Then the fundamental solution \1/ of the partial differential equation 
E^(zi, z n ) = S(zi, z n ) with S = L p g satisfies the identity: 

(5) = (r*tf) * (TJtf). 

We seek the real fundamental solution ^ = \J/ P)9 , since the hyperbolic oper- 
ator L P) g is real. Using the Fourier transform with the generator i by real 
variables z±, z n we infer that 

(6) F(* M )(x) = [F(T*% >g )(x)}[F(Tt% >g )(x)} 

= [-i(EU Xji 2j +Y;] =p+1 x j i :y ii)F(* Pi9 )(a;)][-i(E?=i %+£? =p+1 wl)^,,)^)], 
since 

F(T*¥ M )(:r) = Ej=i F{99 M /dz i )i 2i + E" =p+ i 
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= -i(E]f=i ffjijy +H]=p + iX j i 2 jk)F(^ P;q )(x) 

and analogously for T*, where x = (xi, x n ), Xx,...,x n G R (see also §2). 
For the function 

(7) P(x) = EU x) - E] = p+i Xj with p > 1 and q > 1 the generalized 
functions (-P(s) + iO) A and (-P(x) — iO) A are defined for any complex number 
A G C = R©iR (see Chapter 3 in j(i ). The function P x has the cone surface 
P(zi, z n ) = of zeros, so that for the correct definition of generalized 
functions corresponding to P x the generalized functions 

(8) (P(x) + ciO) A = lim < cei e^o(^(a;) 2 + e 2 ) A / 2 exp(iA arg(P(x) + ice)) 
with either c = — 1 or c = 1 were introduced. Therefore, the identity 

(9) F(%, q )(x) = -iXU x ) ~ E] =P+l x])[F(^ Ptq )(x)r or 

(10) = -l/{P{x) + ciO) follows, where c = -1 or c = I. 

The inverse Fourier transform in the class of the generalized functions is: 

(11) F- x ((P(x) + ci0) x )(z 1 ,...,z n ) = exp(-7rcgi/2)2 2A+n 2 2A+n 7r™/ 2 r(A + 
n/2)(Q(z 1 ,...,z n )-ciO)- x -^/[T(-X)\D\^} 

for each A G C and n > 3 (see §IV.2.6 [5]), where D = det(gj^) denotes 
a discriminant of the quadratic form P(x) = y^"^_^ Qj^.XjXk, while = 
Z)?jfc=i g-'^XjXk is the dual quadratic form so that Y%=i 9^ k 9k,i — $i f° r & U 
j, l; 5j = 1 for j = / and 5/ = for j ^ I. In the particular case of n = 2 the 
inverse Fourier transform is given by the formula: 

(12) F- 1 ((P(x)+ciO)- 1 )(z u z 2 ) = -A- l \D\- 1 ' 2 exp(-7rcgi/2) ln(Q(zi, z n )- 
ciO). 

Making the inverse Fourier transform F _1 of the function — l/(P{x) + iO) 
in this particular case of A = —1 we get two complex conjugated fundamental 
solutions 

(13) ^ m ( Zl , ...,z n ) = - exp(7rcgi/2)r((n/2)-l)(P(^, zJ+dO) 1 -™ 2 '*/ (in 
for 3 < n and 1 < p and 1 < q with n = p + q, while 

(14) ty hl (zi,z 2 ) = 4" 1 exp(7rcgi/2)ln(P(,2i,2:2) + ciO) for n = 2, where 
either c = 1 or c = - 1 . 

Another approach consists in using the anti-derivative operators. The 
hyperbolic operator L P)Q and the delta-function are invariant under the Lie 
group Pi9 (R) or all linear transforms of the Euclidean space R n , n = p + q, 
preserving the scalar product (x,y) Pjq = Y^j=\XjVj — Y^jtp+iXjVj invariant. 
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Thus ty Piq can be written as a composition £(P(x)) of some function £(y) and 
of P(x) given by Formula (7). Therefore, we take the corresponding branch 
of the anti-derivative in the form X^Xy^ = £(P(x)). Applying the Fourier 
transform with the generator i we infer that 

(15) F(%, g )(x) = F(X r X Tl 5) = (P(x) + ciO)- 1 F(5) = (P(x) + dO)" 1 . 
Applying the inverse Fourier transform to both sides of Equation (15) one 
gets Formulas (13, 14). 

Thus the results of §§2-25 over the Cayley-Dickson algebra A v lead to 
the fundamental solution of the hyperbolic operator L p q . This means that 
the approach of §§2-25 over the Cayley-Dickson algebras leads to the effec- 
tive solution of any hyperbolic partial differential equation with constant 
coefficients. Thus Formulas of §§2, 8 with the known function \1/ = \l/ n from 
Formulas 28(10, 11) and (13, 14) of this section give the fundamental solution 
of any first and second order linear partial differential equation with variable 
2-differentiable ^4„- valued coefficients, z G U C A v . 

30. Example. The heat kernel. Each function of the type f(z) = 
P n {z) exp(— t\z\ 2 ) with a marked positive parameter can be written in the z- 
representation due to Formulas 1(1 — 3), where P n (z) denotes the polynomial 
by z of degree n. Therefore, f(z) in the ^-representation is z-differentiable, 
consequently, infinite z-differentiable (see [23l 22j) and 

lim f( m \z).(h 1 ,...,h m )(l + \z\ k ) = 

\z\— >oo 

for each < m, k G Z and Cayley-Dickson numbers hi, h m G A v . There- 
fore, the space E of infinite z-differentiable tending to zero at infinity func- 
tions together with their derivatives multiplied on the weight factor (1 + \z\ k ) 
is infinite-dimensional. Thus it is worthwhile to consider the topologically 
adjoint space E' q of R-linear ^-additive continuous ^-valued functionals 
on E. Elements of E' q are also called the generalized functions. A function 
or a generalized function is called finite if its support is a bounded set. 

The heat partial differential equation reads as 

(l)dv(z)/dz = a 2 Av(z) + f(z), 
where z = z i + ... + z m i m , z , z m G R, 1 < m < 2 V — 1, 2 < v, where 
a > 0, f(z) is a real- valued generalized finite function so that f(z) is zero for 
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zq < (see §16 |37j). We shall seek the generalized solution £ of this equation 
using the technique given above. The generalized function v = £ * f is the 
solution of (1), where 

(2) d£(z)/dz - a 2 A£{z) = 5(z), 

(3) (£*f)(x)= / £(x - z)f(z)dz ....dz m . 

JO JR m 

As usually 5 denotes the 5 generalized function so that 

rxo r 

(4) (8*f)(x)= / 5(x-z)f(z)dz ....dz m = f(x) 

JO JR m 

for each continuous bounded function /. If / is (super-) different iable and 
bounded in each domain {z : < zq < T} for < T < oo, f(z) = 
for z < 0, then the solution v is also (super-) different iable in the domain 
z > as it will be seen from the formulas given below. Let us seek the 
generalized solution £ in the form £(z) = w(zo)e u ^ z \ where w and u are 
unknown real-valued functions to be determined. Calculating derivatives of 
£ and substituting into Equation (2) one gets: 

(5) e u ^{w'(z )+w{z )du(z)/dz }-a 2 e u ^w(z ) £™ 1 [{du{z)/dz J ) 2 +d 2 u{z)/dzj] = 
5(z), consequently, 

(6) (l/tw(zo)X(zo) = -du(z)/dz +a 2 j:T =1 [(du(z)/dz 1 ) 2 +d 2 u(z)/dz 2 } + 
e- u ^(l/w(z ))5(z). 

Take now any sequence of continuous non-negative functions f] n with compact 
supports U n such that U n+ i C U n for each n, with f] n U n = {0}, 

(7) / r] n (z)dz ...dz m = 1 

for all n, tending to 5 on the space of continuous functions p(z) on Ri © 
... © Ri m with the converging integral J Rm +i \p(z)\ 2 dz ...dz m < oo: 
(8) lim.yj—j.Qo / R m+i p(z)rj n (z)dz ...dz m =p(0). 

Therefore, we get that on Ri © ... © Ri m \ {0} for z > the following 
equation 

m 

(9) (l/w(z ))w'(z ) = -du(z)/dz + ^^[(du^/dzj) 2 + d 2 u(z)/dz 2 } 

need to be satisfied. The left side of (9) is independent of z — z , hence the 
right side is also independent of z — z . The partial differential operator 
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{du(z)/dz +a 2 J2 1 jLi[(du/dz j ) 2 +d 2 u(z)/dZj]} acting on u is of the second 
order. For each Cayley-Dickson number z G A r the identity z 2 = 2z Re(z) — 
\z\ 2 is satisfied, particularly, M 2 = — \M\ 2 for each purely imaginary number 
M G A r . Therefore, a function u may be only a polynomial by real variables 
Zi,...,z m of degree not higher than two. On the other hand, the Laplace 
operator A and the 5 function are invariant relative to all elements C of the 
orthogonal group m (R) acting on variables zi, z m . Each m (R) invariant 
real polynomial P of the second order has the form a{z\ + ... + z 2 n ) +/3, where 
a and (3 are two constants independent of zi, z m . Thus u as the polynomial 
of zi, ...,z m may depend on \z — z \ 2 — z\ + ... + z^ only. The latter sum 
of squares can be written in the ^-representation with the help of Formulas 
1(1 — 3). This means that 8 has the form: 

(10) E{z) = w(zo) exp{a(zo)(zf + ... + z^) + P( z o)} and Equation (9) 
simplifies: 

(11) (l/w(z ))w'(z ) = -[da(z )/dz }(z 2 + ... + 4) - [d/3(z )/dz ] 

m 

+a 2 a(z ){2m + a(z ) ^j}- 

j'=i 

We can denote w(z )e^ z °^ by w(z ) again and take without loss of generality 
that /3 — 0. The left side of (11) is independent of zi,...,z m , hence terms 
with \z — zq\ 2 in (11) are canceling: a~ 2 (z )[da(z )/dzo\ = 4a 2 . The latter 
differential equation gives a(zo) = — l/(co + 4a 2 ^ ), where Co is the real 
constant. Substituting this a into (11) one gets: 

(12) (l/w(z ))w'(z ) = a 2 a(z )2m. 

Together with Condition (2) this gives Cq = and the heat kernel S: 

(13) £{z) = ^(^o)[2a(7r^ ) 1 /2]-™exp{-|^ - z \ 2 / '[4a 2 z ]} 
and the solution 

(14) v = E*f, 

where 6{zq) = 1 for zq > and 9(z ) = for z < 0. 

If use anti-derivation operators the solution has the form 26(6 — 8) suppos- 
ing that a solution S is real-valued on real- valued test functions ou, [£, u>) G R, 
where Ti = T, TT = — a 2 A rn (see §) and T 2 = d/dz . Therefore, 

(14) a 2 S = -l r l r (dS/dz - 5). 
Making the Fourier transform F = F ZUmmmiZm by the variables zi,...,z m with 
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the generator i of both sides of Equation (14) one gets for suitable branches 
of the anti-derivatives 

(15) a 2 F(£)(z^x ll ... 1 x m ) = [a 2 Y™=ix]Y\dF(£)/dzu-5(z )). 
Solving the latter ordinary differential equation one finds F(£)(zq, x±, x m ) 
and making the inverse Fourier transform by the variables x±, x m G R one 
gets Formula (13). 

31. Example. Wave operator. In this section the fundamental solu- 
tion £ = £ n of the wave operator is considered: 

(1) U£(t,x) = S(t,x), where 

(2) 11/ = d 2 f/dt 2 — Af denotes the wave (d'Alambert) operator with 

n 

(3) Af = J2d 2 f/dxl 

3=1 

where t > 0, xi, x n G R. We make the change of variables putting t — Z2, 
Xj = z 2j+2 for each j = l,...,n, z = z i + z 2 i 2 + ■■■ + z 2 v_ 2 i 2 v_ 2 G Ai, v , 
z ,...,z 2 v-i G R, r = 1. We consider the case n = 3 and v = 4 so that 
.4.1,4 is isomorphic with the octonion algebra A3 = O. Let us seek £ in 
the class of the generalized functions in the form £(z) = Q(z2)f{z), where 9 
and / are some generalized functions to be calculated, / may depend only 
on -22,-24, ...,z 2n +2- D'Alambert's operator JJ is invariant relative to any R- 
linear transformations A from the Lie group Oi jTV (R). Elements of the group 
Oi jn (R) are characterized by the condition A l GA = G, where G denotes 
the square (n + 1) x (n + 1) diagonal matrix G = diag(l, — 1, —1), the 
transposed matrix A is denoted by A 1 . This means that the wave operator 
TJ is invariant under change of variables £ = (-22,-24, ...,z 2n+2 )A for any A G 
Oi j7l (R). Making the differentiation of £ one gets the differential equation: 

(4) JJ£ = (d 2 6/dzl)f + 2(d9/dz 2 )(df/dz 2 ) + em = S(z). 

The 5-function 5 (z) is also invariant relative to all transformations of 
the Lie group Oi >n (R), since Sg = g(0) for each continuous function g with 
J Rn +i \g(z)\ 2 dz 2 ...dz 2n+2 < 00. On the other hand, Equation (4) implies that 

(5) d 2 6/dz 2 = -[2{d9/dz 2 ){df/dz 2 ) + eUf - 5{z)\/f{z) 

for each z G A r , v , when f(z) 7^ 0. The left side of Equation (5) may depend 
only on z 2 , consequently, the right side of (5) is independent of -2 4 , ...,-2 2n+ 2- 
In view of Formulas 29(2, 3) with p — 1 and q = n we get the operators a — T 
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and o\ = Ti with Q = up to the enumeration of the variables. Therefore, 
one gets the functions \I/i jn (see Formulas 29(13, 14)) over the Cayley-Dickson 
algebra A v . But due to the Oi jn (R) invariance of the generalized function £ 
we infer that it is necessary to take the Oj )n (R) invariant polynomial P(y) = 
{Vl-Y. U y! j+2 )- Thus we put £ = 6(z 2 )f(z) with f(z) = u(z 2 -J2] =1 4 +2 ), 
where u is some generalized function. Substituting u instead of / into (5) 
one gets the simplified differential equation. If suppose that dO/dz 2 = for 
z 2 > 0, then d 2 9 jdz\ = and Equation (5) leads to the differential equation 

(7) 4w( 2 ).(l, l)(r})r] - 4u'.l(<n) = S(z)/c, 

where r\ = r](z) = z\ — z\ — z\ — z%, 0(z 2 ) = c = const for z 2 > 0. Choose any 
sequence of z-differentiable functions g n (z) with compact supports converg- 
ing to the 5-function as in §§24 and 30, when n tends to the infinity. Each 
function g n (z)/r](z) has poles of the first order at points z 2 = \z\ + z 2 + z 2 ] 1 / 2 
and z 2 = —[zl + zl + z^] 1 ^ 2 . Making the substitution p = u'.l in (7) and For- 
mula 3(10) [20j with the right side Q(z) = g n (z)/r)(z) and b(z) = —l/r)(z) we 
obtain the integral expression for the solution p n of the differential equation 

(8) P'n-H 7 !) ~Pn(v)/v = 9n{z)/(Acri). 

To evaluate the appearing integrals it is possible to use Jordan's Lemma (see 
§2 in [24]) over the octonion algebra isomorphic with A\^. The evaluation 
of the integrals (see §3 also) with the given functions can be reduced to the 
complex case, when a and x belong to the same complex plane Cm- Applying 
Jordan's lemma one deduces the expression for p n and the limit function 
p(i]) = 5'(f])/(2^c) + K, where K is a constant, since rj £ R. Therefore, 
u(rj) = J p(i])dr]. Thus one infers the fundamental solution 

(9) £ 3 (z) = 6{t)5{t 2 - M 2 )/(27r) 

and the generalized solution £ 3 * s of the wave equation 

TJ/ = s, where s = s(z) is & generalized function or particularly a z- 
differentiable function. The delta generalized function 5(P) of the quadratic 
form P(x) = x\ + ... + x 2 — x 2 +1 — ... — x 2 +q is described in details in §IV.2 
®. 

32. Helmholtz' operator. 

When (3 ^ with Re{t*f3) = for each j (see §8), for example, = 
Po + ikfii with real (3q and Pi and k > n, then 
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(1) A = A„ + \j3\ 2 is Helmholtz' operator. The corresponding operator 
T is given by Formula 28(2). 

For an arbitrary real non-degenerate quadratic form P(x) generalized 
functions (c 2 + P + biO) x with c > 0, b = 1 or b = -1, A £ C = R ® Ri, are 
defined as: 

(2) (c 2 + P + M0) A = lim 0<e _K)(c 2 + P + 6dPi)\ 
where Pi is a positive definite quadratic form. 

Some special functions are useful for such equations. Bessel's functions 
are solutions of the differential equation 

(SI) z 2 d 2 u/dz 2 + zdu/dz + {z 2 - X 2 )u = 0, 
where A and z are complex. Bessel's function of the first kind is given by the 
series: 

oo 

(S2) J x (z) = J2 (-l) m (z/2) 2m+x /[m\T(X + m + l)\, 

m=0 

where z and A £ Ci. Then 

(S3)/ A (^)=exp(-7rAi/2)J A (k) 
is called Bessel's function of the imaginary argument. Other needed functions 
for non- integer A are: 

(SA) N x (z) = [J x (z) cos(ttA) - J_ A (z)]/sin(7rz), 

(55) HW{z) = J x (z) + iN x (z), 

(56) H( 2 \z) = J x (z)-iN x (z), 

(57) K x {z) = 7r[/_ A (z)-/ A (z)]/[2sm(7rA)] = vr exp(i7r(A + l)/2)^ 1) (iz)/2 
with the complex variable z and non-integer complex parameter A. For inte- 
ger A values of these functions (5*4 — 7) are defined as limits by A £ Ci \ Z. 
The functions and are also solutions of Bessel's differential equation 
(5*1) and they are called Hankel's functions of the first and the second kind 
respectively, K x (z) is known as Mcdonald's function. The functions I x {z) 
and K x (z) are linearly independent solutions of the differential equation: 

{S8) z 2 d 2 u/dz 2 + zdu/dz - {z 2 + X 2 )u = 
(see about special functions in [T81 1ST]). 

Analogously to §28 using Formulas 19.4(C4— CI) and 14(3, 4) or Theorem 
23 for a fundamental solution \l/ n of the equation 

(3) A^ n = 5, 



70 



where A is Helmholtz' operator, we get the identity 

(4) F(¥ B )(x) = [c 2 - (E] =1 x 2 )][F^ n )(x)} 2 or 

(5) F(* n )(x) = l/[c 2 - (E"=i^ 2 ) + MO], 

where c > 0, c = The Fourier transform of the generalized function 
(c 2 + P(x) + 6i0) A by the real variables x = (xi, ...,x n ) with the generator i 
is: 

(6) F[(c 2 + P(x) + bi0) x }(y) = 
2 A + i (27r) n/ 2c A + (n/ 2 )^ A+W2)[c(g(y) _ biO)} /[T(-X)D l/ 2 (Q(y) - biO)W 2+n /% 

where D = det(gj t k) denotes a discriminant of the quadratic form P(x) = 
J21 t k=i93,h x j x k, Q(y) = Y^jk=i9^ kx 3 x k is t ne dual quadratic form so that 
El=i9 j ' k 9k,i = Sj, 5j = 1 for j = I and 5j = for j ^ I (see §IV.8.2 [6]). 
Mention that D 1 / 2 = |D| 1//2 exp(— qni/2) if the canonical representation of 
the quadratic form P has q negative terms. 
Another formula is: 

(7) F[(c 2 +P(x) + bi0) x }(y) = 
2^H270 n/2 exp(-&g7ri/2)c A +^ 
Hb7ri/2)H^}\ n/2) [c(Q^ 
where = 1, = 2, 

(P$,<j>) = ! P>0 P x (j)dx l ...dx w 

(P x ,(p) = j P<0 P x (pdx 1 ...dx n . 
The functions (P + 6i0) A and P x and P x by the complex variable A are 
regularized as it is described in [B] with the help of their Laurent series in 
neighborhoods of singular isolated points A such that after the regularization 
only the regular part of the Laurent series remains. The functions (P + bi0) x 
with b = 1 or b = — 1 have only simple poles at the points A = — n/2, 
-(n/2)-l,...,— (n/2) — k,..., where k = 1,2,... is a natural number. Using 
formula (6) with A = — 1 and P(x) = — (x 2 + ... + x 2 n ) one gets the fun- 
damental solution \I/ n , where (F~ 1 g)(x) = (2n)~ n (Fg)(—x). Particularly, 
^3(2;) = — exp(6ci|x|)/(47r|x|), ^2( x ) = ~^Hq\c\x\)/A or its complex conju- 
gate iHq 2 \c\x\) /4, where Hq denotes Hankel's function, j = 1,2. 

33. Klein-Gordon's operator. 

Consider (3 and tj as in §8 with Re(t*(3) = for each j, c — |/3| > 0. Take 
the operator 
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(1) A = L M + c 2 , 

where L pq is the hyperbolic operator as in §29. For p = 1 and q = 3 the 
operator Aq is called Klein-Gordon's operator. From Formulas 32 (C4 — C7) 
and 14(3,4) or Theorem 23 we infer that 

(2) F(? n )(x) = [c 2 - (EU x j ~ E" =p+1 x 2 )][F(vI/ n )(x)] 2 or 

(3) F(* n )(x) = l/[c 2 - (£f =1 xf - E? =p+ ix 2 ) + MO]. 

Then Formulas 19.4(6) or 32(7) with A = — 1 and P(x) = —x\ — ... — x p + 
x p+1 + ... + x 2 , n = p + q, give the fundamental solution $ M of the equation 

(4) (L M + c 2 )^ = 5, 

where (F^ 1 g)(x) = (2ir)~ n (Fg)(—x). There are two R-linearly independent 
fundamental solutions, so their R-linear combination with real coefficients 
ol\ and a 2 such that a.\ + a.2 = 1 is also a fundamental solution. 

34. Remark. Certainly, more general partial differential equations as 
30(1), but with d l v/dz l , I > 2, instead of dv/dzo can be considered. It 
is worth to mention, that alternative deductions using Formulas 7(1) and 
27(11,14) can be used instead of 8(1) and 19.4(C1 - C7) in §§30 and 31 
providing also u(z) = a(zf + ... + z^) + (3 and f(z) = u{z\ — YTj=\ z h+2) 
with the help of the symmetry Lie groups O m (R) and Oi n (R). Indeed, 
Functions P(x) x satisfy Condition 7(P3) for any real A, where P(z) = z\ + 
... + z 2 p - z 2 p+l - ... - zl, 1 < p < n < 2 V - 1 - 1, 1 < v, since (dP(z) x /dz).h = 
P (z) x ~ l \2Re{j](z)h) , where z G A v , r](z) = z\%\ + ... + z n i n for p = n, while 
r](z) = Zii 2 + ... + z p i 2p + z p+1 ili 2 ( p+ i) + ... + z n i\i 2n for p < n. The function 
r](z) can be written in the z- representation due to Formulas 1(1 — 3). 

Formally the case of the hyperbolic operators L VA + c 2 and their funda- 
mental solutions is obtained from the elliptic operators A n + c 2 with c > 
by the change of variables (xi, x n ) >-)■ (xi, ...,x p ,x p+ ii, ...,x p+g i), where 
n = p + q, Xj G R and xf\ G C = R © Ri for each j, since the quadratic 
forms P may be with complex coefficients and their Fourier transforms can 
be considered as in [B|. At the same operators o or T for these particular 
operators L p q + c 2 and A n + c 2 can be written over the complexified algebra 
(A r )c instead of the Cayley-Dickson algebra A v , 2 < r < v (see §2 above). 
For this we take in Formula 2(8) i instead of w* so that 

(1) *f(z) = E P j= i( d f( z )/ dz ^J + EJ =p+1 (9/(z)/%)ii*, consequently, 
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(2) ( c + a)(c-a)f = (L M + c 2 )f and 

(3) (c + a)(c + a)*/ = (c + a)*(c + a)f = (A n + c 2 )f. 
Let 

(4) H C)P)g = L m + c 2 , where c G R, L„ i0 = A n , 1 < p < n, q = n - p, and 
let ^ c ,p,q be a fundamental solution of the equation 

(5) ^c,p,q^c,p,q — 5. 

Then due to Identities (2, 3) a fundamental solution £ = £p+ a of the equation 

(6) (er + (3)£ = 5 can be written in the form: 

(7) £ p+ff = (a + /3)** c , n ,o, where (3 e A r , \(3\ = c, Re(t*f3) = for each j, 
tj = i* for 1 < j < p, tj = ii* for p < j < n. Moreover, 

(8) £ c+a = (C- (?)^c,p,q- 

Therefore, we infer a solution of the equation 

(9) (tr + P)f = g in V(R m ,A r ) or in the space V(R m ,A r )*: 

(10) / = £/3+a * g. From (2, 3) we deduce a fundamental solution V of the 
equation 

(10) A V = (a + /3)(o"i + /3 X )V = 5'm the convolution form: 

(11) V = £ a+ p * £ ai +h, 
since 

(12) A V = (a + 0)(a, + A) (((a + /3)** c , n , ) * ((^ + /?i)** c ,„, ) 

= (((a + f3)(a + Pi)*Vc, n ,o) * (((cTi + + 0i)**c,n,o)) = 5*5 = 5. 

Particularly, 

^c,p, q = -£ c +a*£-c+a = ((c- o)^ c , n ,o) * ((c+ er)*^) , that can be lightly- 
verified after the Fourier transform by real variables with the generator i, 
since by Formulas (1, 2) the operator a and its anti-derivative operator X CT 
correspond to the signature (p,q) and F(aa*g) = — \z\ 2 F(g) for any g G 
V{W\A v )l 

Knowing a fundamental solution it is possible to consider then a boundary 
problem (see also [T3| 137]). 

35. Partial differential equations of order higher than two. 

The fundamental solution \l/x m +/3 of the equation 

(1) 2 2m ,/3\l/r m +/3 = 5, 
where 

(2) H 2m ^ = (T m + /3)(T m + /3)* 
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can be obtained using decompositions with the help of operators of the first 
order T fc + f3 k by induction, if such decomposition exists (see §§10-14 above). 
Suppose that this decomposition is constructed 

(3) (T m + 0)f(z) = (T m + /3 m )[...[(T 2 + (3 2 )[(T m + (3 m )f(z)}]...], then the 
fundamental solution can be written as the iterated convolution 

(4) ^ Tm+/3 = [...[[(T m +/3 m )^ m ]*[(T m _ 1 +/3 m _ 1 )^ m _ 1 ]]*...]*[(T 1 +/3 1 )^ 1 ], 
where Sj denotes the fundamental solution of the elliptic second order partial 
differential equation 

(5) AjSj = 5, with 

(6) A j = (r j + p j )(r j + p j y. 

The fundamental solutions were written above in §§2 and 28. Indeed, 
using Equalities 4(7 — 9) by induction we have 

(7) T, s (---(am tkm al l _ ljkm _ 1 )al l ^ 2jkm _ 2 )...)al tk Ja lt iJa 2 ,i 2 )...)a mj i m 
= Re(a* mkm a mi i m )...Re(a* 1M ai :h ), 

where J2 S denotes the sum by all pairwise transpositions (ki,li),...,(k m ,l m ), 
dj tk G A v . Therefore, 

(8) S 2m , = 

[-[[(T m + f3 m )(T m + f3 m yS m ] * [(T m _! + /3 m _i)(T m _i + f3 m ^yS m ^]] * ...] * 
[(Tx + /?i)(Ti + = [...[<J * 5] * ...] * 5 = 5. 

Vice versa if the fundamental solution $x m +^ is known, then we get the 
fundamental solution £™ of the equation 

(9) (T m + p)S$ = 5 as 

(10) Sf = (T m + /3)*^ T m +/3 in accordance with (2,7). Moreover, the 
equation 

(11) A m+k f = g with A m+k = (T™ + ft)(T* + &) 

in X>(R n ,^) or in the space P(R n , A„)*, where n is a number of real vari- 
ables, 2 <n <2 V , has the fundamental solution V m+k : 

(12) V m+fc = ^Tf+^i * ^t*+/3 2! where 

(13) frf+ySi = (Y™ + ^i)*^ r Tf+ / 8i 

denotes the fundamental solution of the equation 

(14) (T™ + /3i)£r™ +/3i = 5, consequently, 

(15) / = V m+k * g is the solution of Equation (11). 

For example, the fourth order partial differential operator 
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AJ(z) = EU ^fW/dz* ~ E"= P+ i d*f(z)/dz* 
can be decomposed as the composition of two operators of the second order 
T 2 and T\ formally as a and o\ in 2(8,9) with the substitution of df/dz 2 rj 
on d 2 f ' jdz\rj so that in accordance with Theorem 10 this operator A 4 can 
be presented in the form given by Formulas (2, 3, 11). 

On the other hand, fundamental solutions of and L k v q and for 
certain other second order partial differential operators are known. So com- 
bining them with operators of the form T™ 1 ...T™ k permits to consider fun- 
damental solutions of many partial differential operators of order higher than 
two as well. 

Thus knowing fundamental solutions of the corresponding first or sec- 
ond order operators permits to write fundamental solutions of higher order 
partial differential operators considered above with the help of the iterated 
convolutions in a definite order prescribed by the induction process. 

36. Non-linear partial differential equations. 

We consider the differential equation 

(1) {T m + P + f(y)T)y = g, 

where T m + f3 is a partial differential operator as in Formula 10(13) of order 
m, f(y) is a A v differentiable function, y = y(z) is an unknown function, 
f(y)Ty := J2]=o[f(y)-(dy(z)/dzj)}4>*(z). Suppose that a fundamental solu- 
tion £r m +/3 of Equation 35(9) for the operator (T m + f3) is known. Find at 
first a fundamental solution y = V of (1) with g = 5. Then 

(2) (T m + (3)V = 5 - /I, 

where fj,(z) = f(y(z))Ty(z). The ant i- derivative gives 

(3) w(y(z)) = T J(f(y)Ty)dz = J f(y(x))dy(x) = J f(y)dy = T / p(z)dz, 
then 

(4) y = w- 1 ( yT J ia][atz] fi( y x)dx), 

where w~ l denotes the inverse function. On the other hand, 

(5) y = £r m +/3 * (5 — = Iy™ +/3 (5 — /i), when T m is either of the first 
order for m — 1 or is expressed as a composition of operators of the first 
order, 

(5.1) T m + p = (Ti + P 1 )...(T m + f3 m ) so that 

(5.2) X T ™ + /3 = X Tm+/3 m...X Tl+/3 i, 
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consequently, (4, 5) imply the equation: 

(6) Sym+p *(6-fi)= Xy m+f3 (5 -n) = w~ 1 {y X 7 «[ [0itz] v{x)dx) or 

(7) w(£ym +l3 *(8-fj,)) = y / 7 «| [a M fj,(x)dx. 
We have the identity 

(d(£ * *)/dz p ), 4>) = -((£**), <90/dz p ) = ((£ * (dV/dzp)), 4>) 
with a generalized function 

Therefore, differentiating (6) by zo,...,z n , we infer that: 

(8) {E"=o[^T'"+/3 * (d/i(z)/^)]</>*WH 

(n + l)-^^^- 1 ^)/^)-^)-^^)/^))]^^^ = T£ Tm+/3 or 

(9) T(X Tm+/3 /i) + 

(» + l)- 1 E-, fc= o[(^- 1 (OM)-(^W-(^W/^))]^(^) = T(X Tm+ ^), 
where 

(9.1) [{{dg{z)/dz)4Az)) ® M*)," 92 ) = W) ® l ^ m ) 
for each real- valued test function u and each j, 8(z) = dg(z)/dz, £ = 
x Lo\, fi(x)dx (see also §§4, 5, 17, 23 and 26). If T and (3 are indepen- 

J \[a,t z ] 

dent of Zj, i.e. = is zero identically on U, then d(Xy m+ ^fj,)/dzj = 

Zy m+/3 (dfi/dzj) (see also Note 23.1). Otherwise the derivative d{Ty m +p\i)jdzj 
is given by Formulas 4(6) and 23(8) and 26(1). The function w is known from 
(3) after the line integration by the variable y, so Equation (8) is linear by 
(d[i(z)/dzj), j = 0, ...,n. It can be solved as in [20]. Calculating \x from (8) 
or (9) we get the fundamental solution: 

(10) V = £ym + p * (5 — n) 

and the (particular) solution of (1) is: 

(11) y = V*g. 

When \£y m+ j 3) uj) is real for each real-valued test function uj or Tf = 
(df /dz )4> (z) with a real-valued function <f>o(z) and the inverse relative to 
the convolution generalized function is known £yL + a such that 

(12) £-jm + p * £or' n +/3 = 5, 
then Equation (8) simplifies: 

(13) T^(^) + (n+l)- 1 E-, fc= o^ +/ 3*[[(^-HO/^-(^(^)-(%(^)/^))]^(^)] = 
J2] =0 [£yL + p * {d£ym + p/dz j ))(f)* j (z), consequently, 

(14) TfM(z) + u(ji) = b(z), 

where := (n+l)- 1 E2fe=o^T-+^*[P^ -1 (O/ rf O-(^(^)-(^i(^)/^fc))]0&(^)] 
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and b(z) = SyL +/3 * (T£ X m +/3 ). 

If equation (1) is solved, then it provides a solution of more general equa- 
tion: 

(15) (r m + + f((r) k - 1 z)(r) k )t = g 

finding £ from the equation {T) k ~ 1 C i = y, where (T) fc denotes the k-th power 
of the operator T. 

If <f>j(z) = ijipj(z) for each j, then functions {vj(z) : j} can be chosen 
real-valued or real-valued generalized functions on real valued test functions 
(see System 4(10, 11) and §23). In such case the equality 

= E",fc=o[(^ _1 (0 / d ■ (dvj (z) j dz k )ij)}i* k ilj k {z)\ 
is satisfied. For Tf = (df /dz )(j) (z) with a real- valued function 4>q(z) these 
equations simplify, since 9.h = /i(z)h for each h G R and z G U and 

[(dw-\t)/dt).iM(z)]. 

Thus the results of this paper over the Cayley-Dickson algebras enrich 
the technique of integration of partial differential equations in comparison 
with the complex field. 

It is planned to present in the next paper solutions of some types of non- 
linear partial differential equations with the help of non-linear mappings and 
non-commutative line integration over the Cayley-Dickson algebras. 
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